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Exact  Detection  Performance  of  Multiple-Pulse  Frequency-Shift 
Signals  in  a  Partially-Correlated  Fading  Medium  with 
Generalized  Noncentral  Chi-Squared  Statistics 

INTRODUCTION 

The  transmission  of  a  succession  of  time-delayed  and/or 
frequency-shifted  signal  pulses  through  a  fading  medium  leads  to 
received  waveforms  which  vary  in  amplitude  and  phase  in  a  random 
fashion  and  with  possibly  complicated  statistical  dependencies. 
The  evaluation  of  the  detection  capability  of  matched  filter 
processing  and  incoherent  combination  in  such  fading  situations 
and  noise  has  been  the  subject  of  many  studies  over  the  years 
[1  -  15].  Some  of  these  efforts  have  been  aimed  at  obtaining 
approximations  to  the  performance  of  the  systems  of  interest, 
while  others  have  yielded  exact  results  in  selected  special 
cases . 

Results  for  partially-correlated  fading  have  been  obtained  in 
[ 3;4 ; 6;7;9; 10; 11 ; 12; 14 ; 15] .  In  particular,  in  [4 ; 6; 10; 11 ; 12 ; 14  ] , 
the  characteristic  functions  of  the  decision  variables  have  been 
obtained  in  closed  form  for  the  case  of  partially  correlated 
fading  of  the  received  signal  pulses  with  Rayleigh  amplitude 
statistics  and  additive  Gaussian  noise.  The  results  in  [12]  are 
approximate;  however,  the  signal  amplitude  fading  is  not  limited 
to  Rayleigh  statistics  there,  but  rather  can  have  a  Chi 
distribution  of  an  arbitrary  number  of  degrees  of  freedom. 

The  case  of  independent  signal  fading  in  noise,  but  allowing 
correlated  clutter,  was  addressed  in  [13].  An  extension  of  [12] 
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to  partially  correlated  signal  fading  in  a  system  with 
normalization  was  solved  in  an  approximate  fashion  in  [14]. 
Finally,  a  case  of  partially  correlated  signal  fading  with 
Rayleigh  amplitude  statistics  in  the  presence  of  noise  and 
K-distributed  clutter  was  solved  exactly  in  [15]. 

Here,  we  will  give  exact  results  for  a  fading  medium  which 
has  both  deterministic  components  as  well  as  random  components 
with  arbitrary  covariance  coefficients  and  additive  Gaussian 
noise.  Also,  frequency  shift  keyed  (FSK)  signals  will  be 
allowed,  with  arbitrary  (non-overlapping)  occupancy  in  time, 
frequency  space.  In  addition,  the  signal  amplitude  statistics 
will  not  be  limited  to  Rayleigh,  but  will  include  the  range  of 
possibilities  inherent  in  the  noncentral  Chi  distribution  of  an 
arbitrary  number  of  degrees  of  freedom.  The  end  result  is  a 
closed  form  for  the  characteristic  function  of  the  decision 
variable  that,  although  complicated  in  appearance,  is  amenable  to 
efficient  computer  evaluation  of  the  detection  characteristics  by 
means  of  a  single  fast  Fourier  transform  ( FFT ) .  No  matrix 
inverses  are  required.  However,  a  novel  expansion  technique  for 
the  detection  probability  is  derived  which  is  efficient  and 
accurate;  a  program  incorporating  this  expansion  is  listed  and 
exercised  for  several  examples. 

Extensions  to  a  more  general  form  of  processor  will  also  be 
solved  but  not  evaluated.  In  particular,  the  characteristic 
function  of  the  most  general  complex  second-order  form  of 
arbitrarily-correlated  complex  Gaussian  random  variables  with 
arbitrary  means  will  be  obtained  in  closed  form. 
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RECEIVER  PROCESSING  FOR  A  SINGLE  PULSE 

In  this  section,  we  will  describe  the  basic  model  of  the 
single-pulse  transmitted  signal,  the  received  signal  with  fading, 
the  received  noise,  and  the  receiver  processing.  The  extension 
to  multiple  signal  pulses  will  be  undertaken  in  the  next  section. 

We  presume  that  the  transmitted  real  signal  s(t)  is 
narrowband  with  low-frequency  complex  envelope  s^(t)  superposed 
on  carrier  frequency  fQ: 

s(t)  =  Re{s^(t)  exp(i2nfQt)}  .  (1) 

The  transmitted  signal  energy  is  then 

E  s  J  dt  [s(t)]2  =  \  J  dt  |s(t)|2  ,  (2) 

where  integrals  without  limits  are  over  (-»,+<»).  The  time- 
bandwidth  product  of  the  single  pulse  complex  envelope  s(t)  is 
arbitrary;  thus,  for  example,  s(t)  could  contain  linear  frequency 
modulation  along  with  rectangular  or  Gaussian  amplitude 
modulation. 

The  received  signal  waveform  is 

Re {r  s(t-td)  exp[i2n(f0+fd)t  +  i6]}  ,  (3) 

where  r  is  a  (dimensionless)  amplitude  scaling,  td  is  a  time 
delay,  fd  is  a  frequency  shift,  and  0  is  a  phase  shift.  Real 
random  unknowns  r  and  0  do  not  vary  with  time  over  the  duration 
of  pulsed  signal  s(t).  Delay  td  and  shift  fd  are  presumed  known 
at  the  receiver.  The  average  received  signal  energy  is,  using 
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(3)  and  (2), 

r2  \  J  dt  |s(t) |2  =  r2  E  .  (4) 

The  received  zero-mean  additive  noise  waveform  n(t)  is 

n(t)  =  Re{n( t)  exp(i2nfQt)}  ,  (5) 

which  is  presumed  stationary  over  the  duration  of  the  signal. 

The  spectrum  of  received  noise  n(t)  in  the  neighborhood  of 
frequency  f  is  flat,  with  one-sided  spectral  level  NQ  watts/Hz. 

The  reference  waveform  employed  at  the  receiver  corresponds 
to  the  matched  filter  to  the  transmitted  signal,  namely 

A  Re{s(t-td)  exp[i2n(f0+fd)t  +  i$]}  ,  (6) 

which  utilizes  knowledge  of  td  and  fd-  The  local  reference  level 
A  and  local  phase  shift  <}»  in  (6)  are  irrelevant  to  the  processing 
employed  here;  that  is,  the  performance  in  terms  of  the  receiver 
operating  characteristics  (detection  probability  PD  versus  false 
alarm  probability  Pp)  is  independent  of  A  and  $,  at  least  for 
this  case  of  one  signal  pulse. 

We  now  define  the  two  analytic  functions 

a(t)  *  r  s(t-td)  exp[i2n(fQ+fd)t+i0]  +  n(t)  exp(i2nfQt)  , 

0{t)  *  A  s(t-td)  exp[i2n(fo+fd)t+i$]  ,  (7) 

which  are  recognized  as  corresponding  to  the  received  waveform 
and  the  local  reference,  respectively.  The  total  output  of  the 
matched-filter  squared-envelope  detector  to  the  received  signal 
and  noise  waveform  is  then  given  by 
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Y  =  |J  dt  a(t)  p*(t] 


r  exp(iG)  A  Jdt  |s(t-td)|2  +  A  Jdt  n(t)  s*(t-td)  exp(-i2nfdt) 


~  2 

2rAE  +  n  +  1  n.  *  p  +  n+  in. 

r  1  r  l 


(8) 


where  we  used  (2)  and  defined  the  zero-mean  processor  output 
complex  noise  variate 

nr  +  i  n^  =  exp(-iO)  A  J  dt  n(t)  s*(t-td)  exp(-i2nfdt)  .  (9) 

The  two  covariances  of  the  received  noise  complex  envelope 
n(t)  are  derived  in  appendix  A;  they  are  given  by  (A-8)  and 
( A- 1 1 )  as 

nft^  n(t2)  =  0  ,  n^)  n*(t2)  -  2  Nq  &(tj-t2)  .  (10) 

Use  of  these  relations  with  (9)  enables  us  to  determine  the  two 
averages 


(nr  +  i  ni)  =  nf  -  ni  +  i  2  nr  ni  =  0  , 


(11) 


nr  +  i  ^l2  =  n2  +  n2  =  2  Nq  A2  J  dt  |s(t-td)|2  =  4  Nq  A2  E  , 


where  we  also  used  (2).  Combining  the  information  in  the  two 
lines  of  (11),  there  follows,  for  the  two  zero-mean  real  variates 
nr  and  n^, 

n2  =  n2  =  2  Nq  A2  E  s  a2  ,  nr  ni  =  0  .  (12) 

These  results  are  independent  of  the  parameters  0,  td,  and  fd. 
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Since  noise  input  n(t)  is  Gaussian  and  operation  (9)  is 
linear,  real  variates  nr  and  n^  are  Gaussian  with  joint 
probability  density  function 


P(nr,ni) 


2nad 


exp 


2  2 
4 +  4 


2a 


(13) 


From  this  result,  we  can  derive  the  conditional  exceedance 

distribution  function  of  matched  filter  output  y  in  (8),  for  a 

given  value  of  random  variable  r.  In  terms  of  the  two  auxiliary 

2 

parameters  p  and  a  defined  in  (8)  and  (12)  respectively,  we 
have,  for  u  >  0,  exceedance  probability 

Prob(y  >  u)  =  Probf(p  +  nr)^  +  n?  >  u]  = 


-II 


2  no' 


J  J 

c 

2n  a 

CO 

n 

[  dr' 

j 

r'  [ 

de 

/u 

-n 

dr'  r' 

exp 

2  2' 

.  (*  -  »)  +  y 


2a 


r'2  -  2nr'cos6'  + 
2a2 


/u 


2  2 
r'  +  u 


2  o' 


(“rl-oM'  <14» 


where  C  is  the  region  exterior  to  a  circle  of  radius  /u  centered 
at  the  origin  of  the  x,y  plane,  and  Q  is  Marcum's  function  [1]. 

The  conditional  probability  density  function  of  y,  for  the 
same  given  value  of  amplitude  scaling  r,  is 

P(u>  =  ”  du  Prob<Y  >  u)  =  ~^2  exp 

2  a 


u  +  u 


2  a 


(15) 
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for  u  >  0.  Therefore,  the  conditional  characteristic  function  of 
the  squared-envelope  detector  output  y,  for  given  r,  is 


00 

f  (l)  =  [  du  exp(i£u)  p(u)  *  — ^  f  du  exp  i£u 

J  2aZ  J„  l 


(16) 


where  we  used  [16;  6.631  4].  The  two  important  parameters  here 
were  defined  in  (8)  and  (12)  according  to 

j  -  2  r  A  E  ,  a2  =  2  NQ  A2  E  .  (17) 

Notice  that,  in  this  c^  e  of  matched  filtering,  the  detailed 
behavior  of  complex  envelope  s(t)  is  not  relevant;  only  the 
total  transmitted  signal  energy  E  is  of  consequence  to  the 
conditional  characteristic  function  (16)  of  output  y  in  (8).  The 
results  in  (16)  and  (17)  are  independent  of  the  particular  value 
of  the  medium  phase  shift  0  encountered  in  (3)  or  (7);  this  is 
due  to  the  processing  method  adopted  at  the  receiver,  namely 
envelope  detection  of  the  matched  filter  output. 

The  cumulants  of  random  variable  y,  for  given  fixed  amplitude 
scaling  r,  are  available  by  developing  the  logarithm  of  (16)  in  a 
power  series  in  i£: 

In  fc(l)  =  -  ln(  1  -  H2a2)  +  Up2  (l  -  H2c2)  1  = 

=  <H)j  -  +  ^2(2a2)3-1)  .  (18) 
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There  follows  immediately  the  j-th  (conditional)  cumulant  of  y  as 

2 

X(j)  =  CM)*  (2a2)j  fl  +  j  -AA  for  j  M  .  (19) 

1  2<j  ' 

In  particular,  the  first  two  cumulants  are 

2 

X ( 1 )  -  2a2  +  p2  ,  X(2)  =  4a4  (l  +  .  (20) 

It  is  convenient  to  define  a  conditional  deflection  criterion 
at  the  squared-envelope  detector  output  y  as  the  ratio  of  the 
difference  of  means,  with  and  without  signal,  to  the  noise-alone 
standard  deviation.  That  is. 


where  we  used  (20)  and  (17).  This  conditional  deflection  depends 
on  amplitude  scaling  value  r,  of  course.  However,  it  is 
independent  of  local  reference  values  A  and  $  in  ( 6 ) ,  as 
expected,  since  performance  measures  should  not  depend  on 
receiver  processor  absolute  levels  or  phase  shifts,  at  least  for 
the  case  of  a  single  pulse. 
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RECEIVER  PROCESSING  FOR  MULTIPLE  SIGNAL  PULSES 

In  this  section  we  will  generalize  to  the  case  where  K  signal 
pulses  are  transmitted,  all  of  which  are  nonoverlapping  in  time, 
frequency  space;  the  results  derived  in  the  previous  section  will 
be  used  freely  in  the  sequel.  The  transmitted  energy  in  the  k-th 
pulse  is  Ek,  where  1  £  k  £  K.  The  k-th  signal  pulse  undergoes 
(dimensionless)  amplitude  scaling  rk  and  phase  shift  0k  in 
transmission  through  the  medium  to  the  receiver.  For  each  k, 
random  variables  rk  and  0k  do  not  vary  with  time  over  the 
duration  of  the  k-th  individual  signal  pulse;  this  is  a  general¬ 
ization  of  (3).  The  k-th  local  reference  waveform  employed  at 
the  receiver  utilizes  constant  amplitude  scaling  Ak  and  phase 
shift  <)»k,  which  can  be  chosen  for  best  performance;  compare  (6). 

It  is  now  necessary  to  generalize  the  definitions  in  (7)  to  K 
pairs  of  analytic  functions,  namely 

ak(t)  =  rk  ^■k(t_td)  exPti2n<  f0+fd)t+iek]  +  -(t)  exP(i2nf0t)  , 

Pk(t)  *  Ak  £k(t-td)  exp[i2n(f0+fd)t+i*k]  ,  (22) 

which  correspond  to  the  k-th  received  waveform  and  local 
reference,  respectively.  A  block  diagram  of  the  signal 
processing  employed  at  the  receiver  is  depicted  in  figure  1.  The 
k-th  signal  complex  envelope  sk( t)  in  (22)  includes  the  time 
delay  tk  and  the  frequency  shift  fk  associated  with  the  k-th 
pulse,  in  accordance  with  the  FSK  pattern  employed  at  the 
transmitter  and  known  to  the  receiver.  The  K  matched  filter 
outputs  are  envelope  detected,  squared,  and  then  summed  to  yield 
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Figure  1.  Block  Diagram  of  Receiver  Processing 


output  y,  which  is  compared  with  a  threshold  for  purposes  of 
declaring  signal  presence  or  absence.  Observe  that  the  matched 
filter  £k(t)  incorporates  the  local  reference  scaling  Ak  and 
phase  shift  $k  and  therefore  accounts  for  (complex)  weighting 
prior  to  the  summation  indicated  in  figure  1.  The  absolute  level 
of  weights  jAk)  does  not  matter  and  can  be  chosen  freely; 
however,  their  values  relative  to  each  other  will  affect  the 
detection  capability  of  the  processor  in  figure  1. 

The  energy  in  the  k-th  (real)  reference  waveform  is,  from 
(22)  and  (2), 

\  J  dt  |pk(t)|2  =  A2  Ek  ,  (23) 

where  Ek  is  the  transmitted  energy  in  the  k-th  signal  pulse.  In 
analogy  to  (17),  we  define  the  K  parameters 

Pk  =  2  rk  Ak  Ek  ,  CT2  =  2  Nq  a2  Ek  for  1  £  k  £  K  .  (24) 
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The  noise  density  level  NQ  is  a  common  constant  applicable  to  all 
K  pulses  because  the  received  noise  is  stationary  over  the  total 
time  extent,  and  it  is  white  over  the  entire  frequency  band  of 
the  K  received  signals. 

The  K  signals  {s^(t)}  are  nonoverlapping  in  time, frequency 
space,  (t,f).  Also,  the  received  noise  is  Gaussian.  Therefore, 
the  K  output  noise  variables  from  the  matched  filters  in  figure  1 
are  statistically  independent  of  each  other;  mathematically,  we 
are  using  orthogonality  relation 

J  dt  !k(t)  st(t)  *  0  for  k  3*  j  .  (25) 

This  is  due  to  nonoverlapping,  in  time  or  frequency  or  both,  of 
all  the  component  FSK  pulses  {s^(t)J,  each  with  their  individual 
delay, shift  tk,fk.  These  observations,  plus  relations  (8)  and 
(16),  allow  us  to  determine  the  conditional  characteristic 
function  of  output  y  from  figure  1  in  the  compact  form 


f 


c 


K  —  1 

(D  =  T~T(l  -  «  2oj;]  exp 

v  2 

_2L,  Pi , 

if  HZ  - - - J 

k=l 1  K) 

k=l  1  -  U  2akJ 

(26) 


It  should  be  noticed  that  this  result  is  independent  of  the 
set  of  phase  shifts  {0^}  encountered  in  transmission;  see  (22) 
and  (24).  The  reason  for  this  development  is  the  method  of 
receiver  processing  adopted  in  figure  1,  namely  envelope 
detection  of  each  of  the  matched  filter  outputs  prior  to 
summation.  However,  (26)  depends  significantly  on  the  amplitude 
scalings  { r^ }  through  parameters  f^k)  in  (24);  explicitly,  the 
alternative  form  of  (26)  is 
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1  -  U  4  N 

o 


exp 


U 


2  2  ~2 
4  rk  Ak  Ek 


1  "  U  4  No  Ak 


E 


kJ 

(27) 


Also,  characteristic  function  (27)  is  independent  of  the  receiver 
set  of  phase  shifts  (♦jc{  used  in  local  references  (22),  again  due 
to  the  method  of  processing  employed  in  figure  1.  However,  (27) 
depends  upon  the  receiver  weights  {A^}  in  a  complicated  fashion. 
In  order  to  determine  a  reasonable  method  of  selecting  this  set 
of  weights,  we  consider  the  following  development. 

The  (conditional)  cumulants  of  system  output  random  variable 
Y,  for  a  given  fixed  set  of  amplitude  scalings  {r^J,  are 
available  by  developing  the  logarithm  of  (26)  or  (27)  in  a  power 
series  in  it.  Reference  to  (18)  and  (19)  immediately  yields  the 
j-th  (conditional)  cumulant  of  y  as 


X  ( j  ) 


for  j  M  . 


(28) 


In  particular,  the  first  two  cumulants  of  y  are 


We  can  now  define  a  (conditional)  deflection  criterion  dR  for 
the  output  of  the  processor  in  figure  1,  exactly  as  was  done  in 
(21)  earlier  for  the  single  pulse  case.  Substitution  of  (29)  and 
( 24 )  then  yields 
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The  absolute  level  of  receiver  weights  {A^}  cancels  out  of  (30); 
however,  the  relative  sizes  of  these  weights  does  affect  the 
value  of  deflection  dR  attainable. 

~2 

If  the  average  power  scaling  of  the  k-th  pulse,  r^,  is  the 
same  for  all  k,  that  is. 


with  equality  only  if  all  the  weights  A^  are  equal.  This  result 

is  not  surprising,  since  there  is  already  a  built-in  dependence 

of  matched  filter  Pk(t)  on  the  transmitted  signal  energy  E^; 

2  ~ 

namely,  (23)  gives  this  energy  as  A^  E^.  Thus,  we  have  the 
conclusion  that  the  "best"  set  of  weights  (A^!,  for  the  case  of 
equal  power  scalings  (31)  in  the  medium,  is  uniform  when  local 
reference  0k(t)  is  chosen  according  to  (22). 

More  generally,  the  average  deflection  follows  from  (30)  as 
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k=l 


r2  A2  E2 
rk  Ak  hk 


N 


K 


k=l 


Ak  sk 


K 

HZ 

k“l 


(33) 


with  equality  if  and  only  if  the  receiver  weights  satisfy 


A 


k 


■  A  r 


2 

k 


h 


for  1  £  k  £  K  . 


(34) 


Scale  factor  A  is  arbitrary,  reflecting  the  fact  that  the 
absolute  level  of  set  {AkJ  is  irrelevant.  It  should  be  noted 
that  the  "best"  receiver  weights  in  (34)  are  independent  of  the 
transmitted  signal  energies  {EkJ.  This  is  again  due  to  the  fact 
that  matched  filter  0k(t)  in  figure  1  and  (22)  already  has  an 
energy  proportional  to  Ek;  see  (23). 

Result  (34)  for  the  weights  (AkJ  requires  knowledge  of  the 

~~2 

average  power  scaling  rk  applied  to  the  k-th  pulse  in 
transmission  through  the  medium.  But  these  average  power 
scalings  may  be  unknown  or  they  may  be  independent  of  k,  the 
particular  pulse  number.  In  this  case,  a  reasonable  choice  of 
receiver  weights  is  simply  to  take  Ak  -  A  for  1  <  k  i  K.  Then 
the  conditional  characteristic  function  in  (27)  of  piocessor 
output  y  reduces  to 


£C(U 


u  4  No  *2  gk 


-1 


exp 


U4A' 


K 

HZ- 

k=l  1 


r2  E2 
rk  Ek 


il  4  Nq  A2  Ek 


(35) 


In  the  (usual)  case  where  all  the  transmitted  signal  energies 
in  the  K  pulses  are  equal,  that  is,  Ek  =  E  for  1  1  k  £  K,  then 
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the  conditional  characteristic  function  in  (35)  reduces  to 


£c<£> 


-  (l  -  1 


H  4  NQ  A2  E 


-K 


exp 


iZ  4  A2  E2 


K 


1  -  i£  4  NQ  A2  E  k=l 


.  (36) 


It  is  now  very  important  to  notice  that  the  only  way  that  the 
amplitude  scaling  factors  (r^l  enter  this  characteristic  function 
is  through  the  single  quantity  (sufficient  statistic) 


K 

S  -d 

k=l 


(37) 


2 

which  is  the  sum  of  all  the  power-fading  variates  r^  on  the  K 
transmitted  pulses.  This  simplified  result  in  (36)  holds  when 
the  following  two  reasonable  conditions  are  satisfied: 


Ek  =  E  ,  Ak  ■  A  for  1  £  k  i  K  .  (38) 

That  is,  the  transmitted  signal  pulse  energies  are  all  equal  and 
the  receiver  scalings  {A^j  in  (22)  are  all  taken  equal.  This 
particular  case  will  receive  most  of  our  attention. 

Since  the  receiver  absolute  scaling  A  in  (38)  is  under  our 
control  and  does  not  affect  detectability  performance,  we  take, 
for  notational  convenience  and  without  loss  of  generality, 

2  1 

PT  =  — - -  .  (39) 

2  N  E 
o 

Then  the  conditional  characteristic  function  in  (36)  of  processor 
output  y  in  figure  1  simplifies  to 

fc(Z)  =  (1  -  i2Z)~*  exp[1  *li2l  &  s]  ,  (40) 
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where  we  used  (37). 

Another  important  observation  must  be  made  at  this  point. 
Suppose  the  characteristic  function  of  random  variable  S,  defined 
by  (37),  is  known;  that  is, 

fs(£)  «  exp(US)  (41) 

is  known.  Then  the  unconditional  characteristic  function  fy(£) 
of  processor  output  y  in  figure  1  is  obtained  by  averaging  (40) 
over  S  and  using  (41);  there  follows 

VE)  -  (‘  -  i2t)"K  fs(r^i2is;)  •  <42> 

This  compact  form  for  the  unconditional  characteristic  function 
of  processor  output  y  depends  critically  upon  being  able  to 
obtain  the  characteristic  function  fg(£)  in  (41)  of  the  power¬ 
fading  summation  S  defined  in  (37).  Therefore,  effort  can  be 
concentrated  on  attempting  to  determine  (41),  either  exactly  or 
by  use  of  several  low-order  moments  of  S;  see  [17]  for  example. 

Of  course,  in  the  process,  the  amount  of  partial  correlation 

.  .  .  2 
between  individual  pairs  of  pulse  power-fading  variates  [rkJ  in 

sum  S  will  come  into  consideration.  In  any  event,  the 

characteristic  function  of  sum  S  in  (37)  is  the  major  item  of 

interest;  given  this  quantity,  the  processor  output  y  in  figure  1 

is  completely  characterized  in  terms  of  its  characteristic 

function  (42),  at  least  when  conditions  (38)  are  satisfied. 

The  more  general  conditional  characteristic  function  in  (27), 

with  arbitrary  |Ek)  and  { Ak } ,  will  be  treated  later,  after  we 

have  introduced  a  detailed  model  of  the  amplitude  fadings  {rkJ. 
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CHARACTERIZATION  OF  FADING 

We  are  interested  in  obtaining  characteristic  function  fs(£) 

2 

of  sum  S  of  power-fading  variates  (rk)  in  (37);  its  use  in  (42) 
will  then  yield  the  desired  characteristic  function  f^(l)  of 
receiver  output  y  in  figure  1.  However,  we  must  first  concen¬ 
trate  on  characterizing  the  fading  which  gives  rise  to  sum  S. 

COVARIANCE  COEFFICIENTS  OF  POWER-FADING  VARIATES  (qk) 

Define  the  k-th  power-fading  variate  qk  as  the  square  of  the 
amplitude-fading  variate  rk  in  (22): 

qk  =  rk  for  U  k  1  K  .  (43) 

The  model  of  fading  that  we  consider  here  is  that  power  scaling 
qk  is  a  sample  of  a  continuous  fading  process  q(t,f),  namely 

qk  =  q(tR,fk)  for  1  S  k  S  K  ,  (44) 

where  two-dimensional  function  of  time  t  and  frequency  f, 

M  2 

q(t,f|-C[c.(t,f)*g„(t,f)|  .  (45) 

Sampling  time  tk  and  frequency  fk  correspond  to  the  time-delay 
location  and  frequency-shift  location,  respectively,  of  the  k-th 
FSK  complex  envelope  signal  £k(t)  employed  in  (22). 

The  number  of  fading  components  in  model  (45)  is  M.  Each 
component  contains  a  deterministic  part  cm(t,f)  and  a  stationary 
zero-mean  Gaussian  part  g  (t,f).  Without  loss  of  generality,  the 
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nonrandom  part  satisfies 

cm(t,f)  i.  0  for  all  m,  t,  f  .  (46) 

Random  parts  { gm ( t , f ) }  of  the  joint  Gaussian  processes  are 
stationary  in  both  t  and  f,  with  covariances 

9m(t,f)  9n(t-r , f-v)  =  Rmn(T,v)  for  1  S  m,  n  S  M  .  (47) 

(It  is  possible  to  generalize  to  nonstationary  Gaussian 
processes  |gm(t,f)J;  however,  each  covariance  Rmn  would  then 
be  a  function  of  the  four  variables  r,v,t,f  instead  of  just 
differences  i,v.)  From  (47),  we  have  property 

R  _(-r,-v)  =  R  ( t , v )  .  (48) 

Fading  model  (44)  -  (45)  does  not  constitute  a  multipath 
medium  but  does  mimic  its  net  effect.  Every  transmitted  pulse 
sk(t)  undergoes  just  one  common  time  delay  td  and  frequency  shift 
f^,  as  indicated  in  (22),  which  are  known  and  utilized  by  the 
receiver.  Rather,  each  signal  pulse  simply  undergoes  a  different 
phase  shift  0^  and  amplitude  scaling  r^  in  (22),  the  latter  of 
which  is  characterized  through  power  scaling  qk  in  (43)  -  (45). 

The  latter  random  variable,  q^,  is  more  general  than  non¬ 
central  Chi-squared  because  the  random  components  {gm(t,f)J  in 
(45)  can  have  unequal  variances  and  can  be  correlated 

with  each  other;  that  is,  we  allow  R _ (t,v)  t  0  for  m  ^  n. 

The  mean  of  power-fading  variate  q^  is  given  by 

%  -  EI  [cm(tk'fk>  +  Rmm(°'°>]  for  1  S  k  S  K  '  (49» 
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where  we  used  (44),  (45),  and  (47).  This  quantity  is  independent 
of  pulse  number  k  if  all  M  deterministic  components  are  zero,  or 
if  they  do  not  depend  on  t  or  f;  that  is,  if  cm(t,f)  is 
independent  of  t  and  f  for  1  £  m  1  M.  Since  is  the 
transmitted  signal  energy  in  the  k-th  pulse,  then 


Ek  c*(tk,fk>  s  DVm  =  deterministic  received  signal  energy 


km 


in  m-th  component  of  k-th  pulse. 


(50) 


M 


M 


E^  >  '  ~  )  Dj^  s  Dk  =  deterministic  received  signal 

m_1  energy  in  k-th  pulse  ,  (51) 


Ek  s  Ejan  =  average  random  received  signal  energy 

in  m-th  component  of  k-th  pulse,  (52) 


M  M 

Er  )  ;  Rnun(0,0)  *  ^  Ekm  s  Ek  "  avera9e  random  received  signal 

energy  in  k-th  pulse  .  (53) 


m=l 


m=l 


The  alternating  (zero-mean)  component  of  power  scaling  qk  is 
then  given  by  (44),  (45),  and  (49)  as 


qk  =  qk  "  qk 


g  [2  c.<Vfk>  W£k>  +  «£<V£k>  -  %,<0.0)] 


(54) 


The  covariance  between  a  pair  of  power- fading  variates  is 


Rkj  =  qk  qj  ' 


(55) 


while  the  covariance  coefficient  between  qk  and  qj  is 
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R 


2si 


PRj  *  TiTsT*  * 
lRkk  RjjJ 


(56) 


Thus,  the  fundamental  calculation  required  for  determination  of 
covariance  coefficient  pkj  is  R^j  as  defined  by  (55)  in 
conjunction  with  (54). 

When  we  substitute  (54)  in  (55)  and  use  the  fact  that 
{gm(t,f)J  are  zero-mean  joint  Gaussian  processes,  the  first-order 
and  third-order  moments  involving  (gm(t,f)}  are  zero,  while  the 
fourth-order  moment  can  be  broken  down  into  a  sum  of  products  of 
second-order  moments,  leading  to  covariance 

■  2mCi  -Owvv  + 

+  ‘WWW  cm<V£k>  cn(tj'fj>  ‘  <57> 

m,n=l  J  J  j  j 

In  particular,  the  variance  of  power-fading  variate  qk  is 


*kk  =  5 1  -  2JZ  + 

m,  n=l 


M 


+  4 


m,  n=l 


Rmn(0'0)  cm(tk/fk)  cn(tk,fk) 


(58) 


The  covariance  coefficient  p^j  between  qk  and  qj  is  obtained  upon 
substitution  of  general  results  (57)  and  (58)  into  (56). 
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A  large  number  of  special  cases  can  be  obtained  from  the 
general  formulation  of  fading  in  (45)  and  (47).  Below,  we  list 
seven  special  cases  that  will  occupy  most  of  our  attention  in  the 
remainder  of  this  report,  and  which  will  consistently  be  referred 
to  according  to  their  number  in  the  sequel. 


SPECIAL  CASE  Is  Zero  means  (cm(t,f)} 


cm(t,f)  =  0  for  1  1  m  1  M  . 


(59) 


M 


R,  .  =  2 


k3  m7n==l  mn  k  j  k  3 


Rkk  -  1C  rL'°'°) 

m,  n=l 


(60) 


SPECIAL  CASE  2:  Uncorrelated  components  fgm(t,f)} 


=  Rrr.«.(T'V)  L  f°r  1  ^  i  M 

mn  mm  mn 


(61) 


M 


Rkj  -  2  i-T  WVVVV  + 

J  m=  l 


M 


+  4 


m=l 


Rmm<Vt1'fk-£j>  cm<V£k>  WV  ' 


'kk 


M  M 

2  EZ  -  4  EZ  <=m(tk.£k) 

m=l  m=l 


(62) 


It  is  very  important  to  note  that  correlated  fading  still  exists 
between  the  received  signal  pulses;  that  is,  covariance  R^j  in 
(62)  is  not  zero,  despite  uncorrelated  property  (61)  in  case  2. 
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SPECIAL  CASE  3 :  Zero  means  and  uncorrelated  components 


pkj 


2  R™(tk'tj'fk'£j)  ' 


skk  *  2  n:  < 

m=  i 


H 


T  ■WWVfj> 

m*8!  J  ■* 


M 


m=l 


mm  ' 


(63) 


SPECIAL  CASE  4 :  Uncorrelated  components  with  identical 
covariances  { Rmn ( t ,  v )  } 

Rmn(T'v)  =  Rn(T'v)  6mn  for  a11  m'  n  *  (64) 

Rkj  -  2M  Rij(tk-tj,£]t-£j)  + 

+  4  j  )  ZZ  cm(tk/fk)  Cm(tj'fj)  ' 

J  J  m- 1  J  J 

o  M  ? 

Rkk  ■  2M  Rn(0,0)  +  4  Rn(0,0)  YU  cm(tk'fk)  *  (65) 
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SPECIAL  CASE  5:  Zero  means  and  uncorrelated  components  with 
identical  covariances 


2M  RhlV-tj.Vfj) 


f 


2M  R^ft^O) 


/ 


tn<Vti>fk~fi> 


lll 


(0,0) 


(66) 


In  this  special  case,  the  covariance  coefficient  between  power¬ 
fading  variates  qk  =  and  Qj  *  Q(hj,fj)  is  the  square  of 

the  covariance  coefficient  between  the  amplitude- fading  variates 
gm(tk,fk)  and  gm(tj,fj);  see  (47).  Due  to  the  identical 
covariances,  (64),  this  latter  covariance  coefficient  is  the  same 
for  every  m  in  the  range  1  £  m  £  M. 


SPECIAL  CASE  6:  Uncorrelated  components  with  proportional 
covariances 

This  is  a  generalization  of  special  case  4;  it  allows  the 
random  components  { gm ( t , f ) |  in  (45)  to  have  different  strengths, 
as  might  be  encountered  in  a  fading  medium.  (Note:  In  the 
following,  the  constant  r^m^  is  the  average  relative  power 
measure  of  the  m-th  random  fading  component;  it  must  not  be 
confused  with  the  random  variable  rk  which  is  the  amplitude¬ 
scaling  on  the  k-th  signal  pulse  in  (22).) 
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IWT''’>  -  RU<t''>>  r(m)  i ^  ,  r*1’  -  1  . 


(67) 


Rkj  *  2  R 


o  M 


M 


.(m) 


+  4  «n«WV£j»  S  r  c”><tR'£k>  cm(tj'£i»  ' 


M 


M 


Rkk  “  2  Rn'°'0)  £11  r(m)  +  4  Rn(°^°)  tZ  r(m)  c^tk^k)  *  (68> 

m»l  m=l 


SPECIAL  CASE  7 :  Zero  means  and  uncorrelated  components  with 
proportional  covariances 


Setting  the  means  {cm(t,f)J  in  (68)  to  zero, 

Rkj  *  2  Ru<Vtj'fk‘fj>  ^  r<m)  ' 


M 


5kk  *  2  R11<0'°) 

m*i 

pkj  " 

[Rn<W£k-£i> 

Rn(0,0) 

(69) 


Notice,  in  this  special  case,  that  are  independent  of  the 

values  of  (r^m^j,  the  relative  power  measures  of  the  M  random 
fading  components. 
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CHARACTERISTIC  FUNCTION  OF  POWER-FADING  VARIATE  q(t,f) 


The  instantaneous  power-fading  process  is,  from  (45), 


M  2 

q(t,f)  *  ZZ  [cm(t'f>  +  gm(t'f)]  “  xT  x  ' 


(70) 


where  Mxi  Gaussian  random  column  vector 

X  -  [cl(t,f)  +  gl(t,£)  .  .  .  cM(t,f)  +  gM(t,f)]T  .  (71) 


We  now  appeal  directly  to  the  very  general  result  derived  in 
appendix  B,  for  the  characteristic  function  of  a  quadratic  form 
and  linear  form,  and  identify  the  quantities  there  according  to 

T 

N  =  M  ,  B  —  I  ,  A  *  0  ,  E  =  [c1(t,f)  .  .  •  CM(t,f)]  , 

M  M 

C  =  Cov(X)  -  gn(t,£)]i  -  [Bmn(0,0)]i  .  (72) 

It  is  important  to  note  that  Mxm  symmetric  covariance  matrix  C  is 
not  a  function  of  t,f,  under  the  stationarity  assumption  (47)  for 
all  M  random  processes  {gm(t,f)J. 

According  to  (B-10),  we  must  solve  the  standard 
characteristic-value  matrix  equation 

C  Q  -  Q  A  ,  (73) 

for  Mxm  eigenvalue  matrix  A  and  normalized  modal  matrix  Q,  where 

T 

A  *  diagjXj  ...  XM]  ,  Q  -  [vx  . . .  VM]  ,  Vm  =  [vml  ...  v^J  , 

(74) 

and  Mxl  column  vector  Vm  is  the  m-th  eigenvector  with  components 
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^ Vmn^ '  1  *  n  *  M*  Then  (B-18)  with  (74)  and  (72)  yields 
deterministic  parameters 


m 


-  V*  E  - 
m 


M 


n=l 


mn 


cn(t,f) 


*■<*'*> 


m 


*  0 


(75) 


for  1  i.  m  £  M.  We  can  now  use  (B-20)  to  obtain  the 
characteristic  function  of  power-fading  variate  q(t,f)  as 


[  m  ,  \] 

— 3 

Int1  -  uwJ 

exp 

it 


M 


m*l 


1  -  i2£X 


m 


(76) 


The  eigenvectors  { Vm }  and  nonzero  means  {cm(t/f)}  enter  this 
result  through  the  terms  { em( t , f ) }  defined  in  (75).  This  general 
result  will  be  simplified,  below,  to  the  seven  special  cases  that 
were  presented  earlier  in  (59)  -  (69). 

By  expanding  the  logarithm  of  general  result  (76)  in  a  power 
series  in  it,  the  cumulants  of  q(t,f)  are  found  to  be 


Xg(P> 


(P-1) 


M 

LZ 

m=l 


m 


P  t 


m 


(t,f)j  XP"1  for  p  M  .  (77) 


For  p  *  1  and  p  *  2,  these  may  be  verified  to  equal  the  results 
in  (49)  and  (58);  the  pertinent  matrix  manipulations  are 
indicated  in  appendix  B  for  the  more  general  case  where  linear 
form  A  *  0.  The  quantities  {XmJ  are  the  eigenvalues  of 
covariance  matrix  C  in  (72),  while  parameters  { em ( t , f ) )  are  given 
by  (75)  in  terms  of  eigenvectors  { Vm }  of  matrix  C  and  mean  vector 
E  *  E(t,f )  in  (72) . 
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SPECIAL  CASE  1:  Zero  means;  see  (59) 


Since  cm(t,f)  =  0,  then  em(t,f)  *  0  from  (75),  and  the 
characteristic  function  in  (76)  reduces  to 


(O  - 


(78) 


SPECIAL  CASE  2:  Uncorrelated  components;  see  (61) 

Here,  from  (72),  we  have 

C  =  diag[Rn(0,0)  .  .  .  R^O^]  .  (79) 

Then,  the  eigenvalue  matrix  A  in  (73)  must  be  identical  to  C,  in 
which  case  (73)  becomes  C  Q  *  Q  C,  where  C  is  diagonal.  This 
forces  Q  to  be  diagonal  also;  finally,  normalization  of  Q  yields 
Q  *  I.  There  follows,  from  (74)  -  (76), 


v 

mn 


m 


*  vt-£> 


and  characteristic  function 


(80) 


VE> 


(ri(l-i2«mm<0'0)) 


-h 


exp 


U 


M 


cm(t'f) 


m*l 


l-i2£R  (0,0) 
mmv 


(81) 


SPECIAL  CASE  3:  Zero  means  and  uncorrelated  components 


Here,  we  simply  set  the  constants  in  (81)  equal  to  zero: 


[  M  / 

\ 

ve>  -  \n( 

1  - i2£ 
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SPECIAL  CASE  4 :  Uncorrelated  components  with  identical 
covariances;  see  (64) 


Now  we  use  (64)  on  result  (81)  to  get 


y«> 


(i-i 


x  -M/2 

i25R11(0,0) j  exp 


f _ U 

ll-i2£R 


11 


(0,0) 


M 


m*l 


(83) 


It  is  important  to  observe  in  this  case  that  only  the  sum  of  the 
squares  of  the  means,  |cm(t,f)},  matters  in  so  far  as  the 
characteristic  function  of  power  scaling  q(t,f)  is  concerned. 


SPECIAL  CASE  5:  Zero  means  and  uncorrelated  components  with 
identical  covariances 


Upon  setting  the  means  in  (83)  equal  to  zero,  there  follows 


(l) 


i  rM/2 

[l  -  i2l  Rn(0,0)j 


(84) 


This  special  case  is  equivalent  to  [12;  ( A— 2 1 ) ]  when  we  take  the 
parameter  m  there  equal  to  M/2. 


28 


TR  10041 


SPECIAL  CASE  6:  Uncorrelated  components  with  proportional 
covariances;  see  (67) 


From  (67),  there  follows 


(0,0)  =  Rn(0,0)  r (m)  6mn  ,  r(1)  =  1 


(85) 


Use  of  this  relation  on  (81)  yields  characteristic  function 


- 


“Tfl-i2^R11 (0,0)r(ra)l  exp  U  JZZ - 

*l'  11  '  m=l  1-i 


c;(^f) 


i2?;R11(0,0)r 


(m)  ‘ 


(86) 


The  cumulants  of  power  scaling  q(t,f)  can  be  obtained  by 
expanding  the  logarithm  of  characteristic  function  (86)  in  a 
power  series  in  i£;  there  follows 


X„(P) 


Art\b/i  t'  M  /m>P  r  *  ♦  ♦  ,  /rnit'  y 

— 3 - ^TT  -  R,i(0»0)  II I  r(  *  +  P  Rn(0,0)  JZZ  r(m)  c“(t,f) 

-1) !  2P  1  11  m=l  11  m=l  m 


(P-1) 


for  pH. 


(87) 


SPECIAL  CASE  7 :  Zero  means  and  uncorrelated  components  with 
proportional  covariances 


Setting  the  means  in  (86)  to  zero,  we  have 


fqU)  -  [rT(l  -  W  Rn(0,0)  r(m)) 


(88) 
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Actual  numerical  examples  displaying  typical  probability 
density  functions  of  power-scaling  random  variable  q,  for  a 
variety  of  different  parameter  values,  will  be  presented  in  a 
later  section. 
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CHARACTERISTIC  FUNCTION  OF  POWER- SUM  VARIATE  S 


The  power-sum  variate  S  is  given  in  (37)  and  (43)  as  the  sum 
of  the  power  scalings  =  q( f ^ )  over  the  total  of  K  pulses 
transmitted  and  received: 


K 

S  =  ZZ  q(tk,fk) 

k=l  K  K 


(89) 


where  we  used  (44)  and  (45).  This  double  sum  can  be  recognized 
as  a  quadratic  form  in  KM  correlated  nonzero-mean  Gaussian  random 
variables.  Therefore,  the  general  approach  in  appendix  B  can  be 
used  directly  to  find  the  characteristic  function  of  S  for  any 
statistical  dependencies  between  the  M  zero-mean  Gaussian 
processes  {gm(t,f)J  and  any  layout  of  the  K  points  {tk,fkl  in  the 
time, frequency  plane.  In  order  to  realize  form  (89),  we  identify 
variate  q  =  S,  constant  N  =  KM,  and  matrices  B  =  I,  A  =  0  in 
appendix  B.  Mean  column  vector  E  in  (B-l)  is  KMxl,  while 
symmetric  covariance  matrix  C  is  KMXKM.  Then  (B-10)  becomes  the 
standard  characteristic-value  matrix  equation  [18;  section  1.13] 


where  KMxKM  modal  matrix  Q  =  [V^  .  .  .  VRMJ  from  (B-ll).  Also, 
there  follows  from  (B-18), 

E  »  QT  E  ,  en  *  vn  E  for  1  S  n  S  KM  ;  aR  =  0  .  (91) 

The  characteristic  function  of  sum  S  in  (89)  then  follows  from 
(B-20)  as 
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(KM  ,  *  KM  E2 

£S<£>  "  [TTl1  -  i2£\,)J  exp  [it  T  _  ?2Unj  .  (92) 

This  general  result  will  now  be  specialized  to  the  seven  cases  of 
particular  interest  here. 

SPECIAL  CASE  1:  Zero  means;  see  (59) 

when  c  (t,f)  =  0  for  1  l  m  £  M,  then  mean  vector  E  is  zero, 
meaning  that  its  components  { en }  in  (91)  are  zero.  The  result  in 
(92)  then  reduces  to 

■  |rr(,  ■ 

Only  the  KM  eigenvalues  of  KMxKM  covariance  matrix  C  of  the  KM 
random  variables  19^^/ )  need  to  be  evaluated  in  this  case. 


SPECIAL  CASE  2:  Uncorrelated  components;  see  (61) 


From  this  point  on,  we  shall  be  interested  in  the  more 
restricted  case  where  component  Gaussian  process  gm(t,f)  is 
uncorrelated  with  (independent  of)  process  g  (t',f')  for  m  /  n, 
regardless  of  the  values  of  t,f  and  t',f';  this  case  was  also 
considered  earlier  in  (61).  This  will  probably  encompass  most 
situations  of  practical  interest;  furthermore,  it  still  allows 
for  correlated  fading  between  the  K  received  signal  pulses.  That 


is,  covariance  in  (62),  between  power 
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need  not  be  zero,  despite  uncorrelated  property  (61)  between 
components  {gm(t,f)}  in  this  special  case  2. 

In  this  case,  a  simpler  and  more  direct  approach  is  possible 
and  is  adopted.  We  begin  by  expressing  (89)  as 

M  (ml 

S  -  ZZ  S(  1  .  (94) 

m=l 

where  we  have  defined  the  M  random  variables 


for  1  1  m  1  M  . 


(95) 


This  superscript  m  notation  is  adopted  in  order  to  readily 
distinguish  the  fading  component  numbers,  1  £  m  £  M,  from  the 
time, frequency  signal  pulse  numbers,  1  £  k  £  K;  see  (45)  versus 
figure  1. 

Now,  it  is  important  to  observe,  for  this  special  case  2, 
that  these  latter  M  random  variables  {S^m^J  are  statistically 
independent  of  each  other,  allowing  us  to  develop  the  character¬ 
istic  function  of  sum  S  in  (94)  in  the  finite  product  form 

f«.(M  -  FT  f(rn)(M  ,  (96) 

b  m=l 

where  the  m-th  characteristic  function  is  given  by  average 

f{m)(Z)  =  exp(us(ra))  ,  (97) 

in  terms  of  quadratic  sum  in  (95). 

In  order  to  ascertain  the  characteristic  function  of  sum 
S(m>,  we  define  Kxi  Gaussian  column  vector  y(m)  according  to 
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,(m) 


[W*!*  +  Wfl>  •  •  •  cm^tK'£K^  +  W^'f-'98) 


Then,  S<m>  in  (95)  can  be  written  in  the  quadratic  form 

s(m)  =  y (m)T  Y(m) 


(99) 


We  now  appeal  to  the  general  results  in  appendix  B  and  identify 
the  quantities  there  according  to  q  =  S^,  N  =  K,  and 


B  =  I  ,  A 


•  0  ,  E  =  E<"»  -  Y<m>  ■=  .  .  .  cm(tK,£K)]T, 


=  p(m), 


C  -  C 


K  K 


(100) 


The  K  diagonal  elements  of  matrix  are  all  equal  to  R  (0,0). 

mm'  ' 

According  to  ( B— 10 ) ,  we  must  now  solve,  for  each  value  of  m 
in  the  range  1  1  m  £  M,  the  Kxk  characteristic-value  equation 


c(m)  Q(m)  .  Q(m)  A(ra) 


(101) 


for  Kxk  eigenvalue  matrix  A*m^  and  corresponding  kxk  normalized 
modal  matrix  Q^,  where 


*  diag|A^m)  . 


X  (n») 
aK 


,  Q<m>  =  [V<m> 


.  V 


(m) 

K 


r(n») 


,(m) 

rkl 


r(m) 

kK 


f. 


(102) 


and  Kxl  vector  is  the  k-th  eigenvector  with  components 

|v<->,.  1  1  j  <l  K.  Then,  ( B— 18)  with  (102)  and  (100)  yields 


(n») 


v(m)T  „(m)  _  (m)  ,t  f 

Vk  E  4—  vkj  cm(tj'  j 


)  ,  aim>  =  0  . 


(103) 
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We  can  now  use  (B-20)  to  obtain  the  characteristic  function 
of  variate  in  (99)  as 


f  ("») 


(O 


exp 


K 


(m) 


« 


k-1  1  -  i2 tk 


(m) 

k 


(104) 


The  desired  characteristic  function  of  power-sum  variate  S  in 
(89)  is  then  given  by  (96),  applied  to  (104): 


*s<£» 


-If 

exp 


M 


(m) 


« 


m=l  k=l  1  -  i2£A 


(m) 

k 


(105) 


This  characteristic  function  of  sum  S  is  a  very  general 
result,  applicable  to  the  case  of  uncorrelated  components 
{ g  ( t , f ) } ;  however,  it  does  require  the  solution  of  M  matrix 
equations  of  the  form  of  (101),  each  matrix  being  of  size  kxr. 
Nevertheless,  this  approach  is  significantly  simpler  than  solving 
the  one  large  KMxrm  matrix  equation  ( 90 ) . 


SPECIAL  CASE  3:  Zero  means  and  uncorrelated  components 


When  the  means  (cm(t,f)}  are  zero,  (103)  yields 
the  characteristic  function  in  (105)  reduces  to  just  the 


=  0  and 
product 


fS<*>  = 


(106) 


35 


TR  10041 


Eigenvalues  A^  ^ are  found  from  Kxk  covariance  matrix 
C(m)  in  (100);  this  solution  must  be  repeated  for  1  £  m  £  M. 


SPECIAL  CASE  4:  Uncorrelated  components  with  identical 
covariances;  see  (64) 


Here,  covariance  matrix  C(m^  in  (100)  is  independent  of  m;  in 
particular,  we  now  have  the  single  Kxk  covariance  matrix 


C  = 


(107) 


This  leads  to  solution  matrices  Q*m)  and  in  (101)  which  are 

also  independent  of  m;  that  is,  (101)  becomes  the  single  Kxk 
characteristic-value  matrix  equation 


C  Q  =  Q  A  . 


(108) 


Thus,  the  eigenvalues  and  eigenvectors  in  (102)  are  independent 
of  m.  However,  the  constants  {e£m)}  in  (103)  still  depend  on  m 
through  their  dependencies  on  means  {cm(t,f)J;  that  is,  from 
(103), 


F(m)  _  VT  (m) 
Ek  vk  11 


m 


(109) 


The  collection  of  all  these  conclusions  enables  us  to  reduce 
the  characteristic  function  (105)  of  S  to  the  compact  closed  form 


fsu> 


-M/2 


exp 


K 

zz 

k=l 


1  -  i2£A. 


(110) 
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where  we  defined  the  K  constants 

M  2 

hk  -  e£m)  for  1  1  k  1  K  .  (Ill) 

m**l 

If  the  number,  M,  of  components  in  fading  model  (45)  happens  to 
be  even,  the  calculation  of  (110)  does  not  involve  a  square  root, 
thereby  further  simplifying  its  numerical  evaluation.  The 
collapsing  of  the  KM  constants  { e ^ }  into  a  smaller  set  of  K 
constants,  by  means  of  the  sums  of  squares  in  (111),  is  the 
analog  of  the  result  in  (83)  for  an  individual  fading  variate 
q(t,f)  in  special  case  4. 

As  a  special  subcase  here,  suppose  that  the  M  deterministic 
components  |cm(t,f)J  in  (45)  are  independent  of  t  and  f;  that  is, 

cm(tk'£k>  “  cm  *or  1  S  k  i  K  .  (112) 

This  means  that  the  constants  are  independent  of  the  locations 
{tk»fkl  of  the  signal  pulses  in  the  t,f  plane,  although  they  can 
still  depend  on  the  component  number  m.  Then,  there  follows  from 
(109)  and  (111), 


(m) 

k 


K 


(113) 


That  is,  the  characteristic  function  of  sum  S  in  (110)  depends  on 
the  constants  { cm }  only  through  their  sum  of  squares.  This 
latter  property  holds  true  only  when  special  subcase  (112)  is 
valid. 
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SPECIAL  CASE  5;  Zero  means  and  uncorrelated  components  with 
identical  covariances 

Upon  setting  the  means  {cm(t,f)J  to  zero,  it  is  seen  from 
(109)  and  (111)  that  characteristic  function  (110)  reduces  to 


(114) 


Now,  only  the  K  eigenvalues  (A^l  of  KxK  covariance  matrix  C  in 
(107)  need  to  be  evaluated.  This  result  in  (114)  is  essentially 
identical  with  [10;  (D-14)].  For  the  special  case  of  M  =  2,  this 
same  result  for  the  characteristic  function  of  S  can  be  shown  to 
follow  from  [11;  (20)]. 


SPECIAL  CASE  6:  Uncorrelated  components  with  proportional 
covariances;  see  (67) 

From  (67)  and  (100),  we  find  that  kxk  covariance  matrix 

C(m)  =  r(m)  c  ,  where  c  .  [RjjtVtj^-fj)]  ,  (115) 

and  r^1)  -  1.  Then  (101)  yields 

Q(ra)  =  Q  ,  A(m)  =  r (m'  A  ,  (116) 

where  Q  and  A  are  the  solutions  to  the  single  kxk  characteristic- 
value  matrix  equation 

C  Q  =  Q  A  .  (117) 
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We  now  refer  to  (103)  to  get  constants 


4m>  =  vk E<m)  -£>kj  wfi 


(118) 


Then  (104)  yields  the  characteristic  function  of  as 


f (m) ( l ) 


“  [FT(i  -  i«V<B’) 


-h 


exp 


K 


(m) 


U 


fel  1  -  i2£Akr(m) 


.  (119) 


Finally,  (96)  gives  the  characteristic  function  of  sum  S  as 


fS<*>  " 


M  K  ,  ,  4.1 

FT  TT  1  -  i2£  Ak  r 

m-1  k*l '  K  ' 


-h 


x  exp 


M  K  e*m) 

- n 57  • 

m*l  k*l  1  -  i2£  A.  r(m,J 


(120) 


Here,  even  if  the  M  deterministic  components  fcm(t,f)}  were 

independent  of  t  and  f,  as  previously  assumed  in  (112),  the 

characteristic  function  of  S  would  not  depend  simply  on  the  sum 
2 

of  squares  £  c  .  To  see  this,  substitute  the  left  member  of 
(113)  in  (120). 


SPECIAL  CASE  7 :  Zero  means  and  uncorrelated  components  with 
proportional  covariances 

Set  the  means  { cm ( t , f ) |  to  zero  in  (118),  at  which  point 
the  characteristic  function  for  S  in  (120)  simplifies  to 
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fs(l)  = 


-  i2l  Xk  r(m)) 


-h 


(121) 


Eigenvalues  {X^}  are  the  diagonal  elements  of  eigenvalue  matrix  A 
in  Kxk  matrix  equation  (117),  where  covariance  matrix  C  is  given 
by  (115). 


CUMULANTS  OF  SUM  S 

For  some  purposes,  such  as  approximating  the  probability 
density  function  or  exceedance  distribution  function  of  S,  the 
cumulants  are  useful;  see,  for  example,  [17].  By  expanding  the 
natural  logarithm  of  the  general  characteristic  function  in  (92) 
in  a  power  series  in  i£,  the  cumulants  of  S  are  readily  found: 

KM 

XS(P)  -  (P-1 )  i  2P_1  ZZ  AP_1(xn  +  pE2)  for  p  M  .  (122) 


The  cumulants  for  the  seven  special  cases  considered  above 
could  also  be  derived.  However,  we  will  only  present  the  results 
for  special  case  6,  namely  uncorrelated  components  with  propor¬ 
tional  covariances;  see  (67).  Again,  expanding  the  logarithm  of 
characteristic  function  (120)  in  a  series  in  i£,  there  follows  a 
slight  simplification  of  (122)  for  the  cumulants  of  S: 


XS(P>  -  (P-1)  I  2p_1  ZZ  r(m)P  ZZ  *£  + 

s  m=l  k=l  K 


i  M  K  , 

♦  pi  2p-‘  ZZ  ZZ 

m-1  k=l  v 


Xk  r 


1P_1 
)  E’ 


(m)lr  *  P(m) 


for  p  2.  1 


(123) 
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CHARACTERISTIC  FUNCTION  OF  PROCESSOR  OUTPUT  y 


With  the  characteristic  function  of  sum  S  in  hand,  we  can  now 
return  to  the  desired  result  (42)  for  the  characteristic  function 
f Y ( £; )  of  processor  output  y,  when  conditions  (38)  are  satisfied. 
Substitution  of  general  result  (92)  into  (42)  yields 


By  expanding  the  logarithm  of  (124)  in  a  power  series  in  i£,  the 
cumulants  Xy(P)  processor  output  y  are  found  to  be  given  by 


(P-1) 


(P) 

!  20- 1 


Mx 

nj 


t1  *  fcj 


P-1 


'n 


for  p  i  1  . 


K(M-2 ) 
(125) 


SPECIAL  CASE  Is  Zero  means;  see  (59) 


Upon  setting  means  {cm(t,f)}  to  cero  in  (124),  the  constants 
{cnJ  defined  in  (91)  become  zero  and  (124)  reduces  to 


f 


Y 


il) 


(1  -  i2MK(‘sM-1) 


1 


n=l 


(126) 


Only  the  eigenvalues  JA.n)  of  KMXKM  covariance  matrix  C  of  the  KM 
random  variables  l9m(  f  j^)  )  need  to  be  evaluated  in  this  case. 
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SPECIAL  CASE  2:  Uncorrelated  components;  see  (61) 


The  comments  in  the  bottom  paragraph  of  page  32  and  the 
sequel  are  relevant  at  this  point  and  should  be  reviewed.  For 
convenience,  we  will  utilize  a  normalized  version  of  covariance 
matrix  C^  defined  in  (100): 


,(m) 


(m) 


R  (0,0) 
mm'  '  ' 


R,„(0,0) 

mm 


K 


(127) 


The  corresponding  eigenvalue  matrix  of  normalized  Kxk 

covariance  matrix  C^m^  is  then  given  by 


(m) 


(m) 


R  ( 0 
mm ' 


(m)  .  .  .  X(m) 

-K 


(128) 


iToT  E  dia9[4 

in  terms  of  eigenvalue  matrix  A^m^  in  (101),  which  now  becomes 

C(m)  Q(m)  *  Q(m>  A(m)  .  (129) 

Normalized  modal  matrix  is  unchanged  from  (101).  The 

relationship  between  the  eigenvalues  in  (102)  and  (128)  is 


\(®®)  =  \(m)  D  (0  0) 

Ak  -k  Rmm(U'U) 

When  we  employ  special  case  2  result  (105)  in  (42),  the 
characteristic  function  of  processor  output  y  is  given  by 


(130) 


fy(U  *  (1  -  i2S)K(lsM~1) 


M  K 

tt  rr 

m=l  k=l 


i  -  m(i  +  x<"") 


-h 


exp 


M  K 


(m) 


u  2E  _ JS _ 

N°  m=1  k=1  1  -  i2s(l  +  ^ 


(131) 
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However,  a  combination  of  (130)  and  (52)  yields 


2E  (m) 

N  Ak 
o 


Rmm(0,0) 
mm'  ' 


v  (m) 
-k 


2Elm  .  (m) 

N  — k 

o 


(132) 


where  we  have  taken  note  that  the  same  signal  energy  E  was 
transmitted  on  all  K  pulses;  see  (38).  That  is,  E^m  is  the 
average  random  received  signal  energy  in  the  m-th  component  of 
any  one  of  the  pulses.  The  quantity  2E^m/NQ  is  a  measure  of  the 
average  random  received  "signal-to-noise  ratio"  in  the  m-th 
component  of  one  pulse.  At  the  same  time,  from  (103), 


2E  (m)  _ 

N  Ek 
o 


K 


(m) 


c_(t 
m' 


j'V  (n 


^ .  \< 


2E] 


(133) 


But  since  cm(t,f)  2  0  without  loss  of  generality  (see  (46)),  we 
have 


(134) 


where  Djm  is  the  deterministic  received  signal  energy  in  the  m-th 
component  of  the  j-th  pulse,  as  defined  in  (50).  Then  (133)  can 
be  expressed  as 


2|  _ (m) 2 

N  Ek 
o 


K 


3  =  1 


,(m) 

kj 


(2DiJ 

lm 

N 

o  J 

E(m) 

-k 


(135) 


The  quantity  2Djm/NQ  is  a  measure  of  the  received  deterministic 
signal-to-noise  ratio  in  the  m-th  component  of  the  j-th  pulse. 
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Upon  combining  these  definitions,  (131)  is  modified  to 


f  (o  =  (i  -  i2E)K<‘>M-1>  frr  rrfi  -  mfi  +  4m» 

Y  lm*l  k*l  1  No  -K 


M  K 


x  exp  ij; 


E(n») 

-k 


:1  k=1  1  -  i2*(l  +  x£m)) 


(136) 


Here,  eigenvalues  Ai  /•••/Ar  are  those  of  the  Kxk  normalized 
covariance  matrix  defined  in  (127).  The  constants  E.£m^  are 

given,  according  to  (135),  by 


F  (m) 

-k 


v(m)  [2Djm 

W  kj  No 


(137) 


The  cumulants  of  processor  output  y  in  this  case  are  given  by 


Xy(P) 
(P-D!  2p‘ 


M  K  (  2Elm  ,.*] 

i  +  -rr2  4  1 

n=l  k=l '  o  k  ' 


M  K 


P  EZ  EZ(l  +  -JT2  ^k”1)  ij'"0  -  K(M-2)  for  p  i  1  .  (138) 

m=l  k»ll  o  K  '  K 


SPECIAL  CASE  3:  Zero  means  and  uncorrelated  components 


When  means  Jc  (t,f)}  are  zero,  then  (50)  and  (137)  yield 
Dk  =  0  and  =  0,  thereby  causing  (136)  to  reduce  to 


M  K 


f  ( o  -  (i  -  i2t>K(*sM_1>  7777  i  -  i2efi  ♦  -jf12  4m))  -<139i 

Y  m«  1  k=  1 1  1  Bo  K  7 
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SPECIAL  CASE  4 :  Uncorrelated  components  with  identical 
covariances;  see  (64) 

Now,  normalized  covariance  matrix  C^  in  (127)  is 
independent  of  m;  in  particular,  we  have  the  single  kxk 
covariance  matrix 


-  Rn(0,0) 


ill(tk~ti/fk~fi  l 


Rn(0,0) 


(140) 


This  leads  to  solution  matrices  and  Av“,,  in  (129)  which  are 

also  independent  of  m;  that  is,  (129)  becomes  the  single  Kxk 
characteristic-value  matrix  equation 


(m) 


C  Q  —  Q  A  ;  A  —  A/R ^ ^(0,0)  diaqfA^, . ■ . , ] 


(141) 


Thus,  the  eigenvalues  in  (128)  and  the  eigenvectors  in  (102)  are 
independent  of  m.  However,  the  constants  in  (137)  still 

depend  on  m  through  their  dependence  on  deterministic  received 
signal  energies  IDjml?  that  is,  from  (137),  now 


F  (m) 

-k 


K 


(142) 


Also,  from  (52),  (38),  and  (64),  we  now  have 

Elm  -  E1  ’W°'°>  *  E  *n(0.0)  =  Eu  ,  (143) 

which  is  the  average  random  received  signal  energy  in  one 
component  of  one  pulse.  (In  this  special  case  4,  we  have  the 
alternative  result  =  E^/M  from  (53).) 
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The  collection  of  these  conclusions  enables  us  to  express  the 
characteristic  function  (136)  of  processor  output  y  in  the  closed 
compact  form 


where  we  defined  the  K  constants 


M  lm\2 

hk  -  EZ  4m) 

K  m-1  K  m 


M  f  K 

rz  fez 

m-l  [j^l 


kj 


(2D. 
jm 


N_ 


for  1  £  k  1  K  .  (145) 


Use  of  (142)  was  made  here.  The  cumulants  of  output  y  are  now 


2E, 


XV(P)  ‘“11  V 

- 1 - — T  =  M  Ol  +  AJ 

-1)!  2P_1  fet'  N0  k' 


(P-D 


K  /  2E. .  i  p~  1 

+  P  HZ  hk  [l  +  IT1  Xkj  -  K ( M-2 )  for  p  i  1 
k*  1  o 


(146) 


As  a  special  subcase,  suppose  that  the  M  deterministic 
components  (cm(t,f)J  in  fading  model  (45)  are  independent  of  t 
and  f;  see  (112).  Then,  from  (50)  and  (38), 


km 


c!(tV,fu)  *  E  C2 


lm 


(147) 


where  D.  is  the  received  deterministic  signal  energy  in  the  m-th 

lift 

component  of  any  one  of  the  pulses.  This  result  allows  the 
constants  in  (145)  to  be  simplified  to 
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-** ' 


(148) 


Thus,  in  this  special  subcase,  the  constants  [h^l  depend  only  on 
the  sum  of  the  received  deterministic  signal  energies  on  all  the 
components . 


SPECIAL  CASE  5:  Zero  means  and  uncorrelated  components  with 
identical  covariances 


Upon  setting  the  means  [c  (t,f)]  to  zero,  it  follows  from 
(50),  (145),  and  (144)  that  the  characteristic  function  of 
processor  output  y  is 


t  (t)  -  (1  -  i2t)K(!,H-1) 


K 

FT 

k=l 


2E, 


11  A 

1  -  i2t[l  +  ~ N~  — k) 


-JsM 


(149) 


Eigenvalues  X^,...,XK  correspond  to  Kxk  normalized  covariance 
matrix  C  given  by  (140),  while  Ej^  is  the  average  random  received 
signal  energy  in  one  component  of  one  pulse.  This  exact  result 
replaces  the  approximation  in  (12;  (8)]. 

In  the  special  case  of  H  =  2,  that  is,  two  components  in 
fading  model  (45),  the  result  in  (149)  reduces  to 


f 


Y 


(l) 


K 


TT 

k=l 


1 


+ 


2E 


11 


N 


(150) 


This  expression  agrees  with  [11;  (24)],  except  for  a  scaling  of 
output  v  by  a  factor  of  2.  There  was  no  need  here  to  evaluate 
multiple  integrals  as  encountered  in  [11;  (21)  -  (22)]. 
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SPECIAL  CASE  6:  Uncorrelated  components  with  proportional 
covariances;  see  (67) 

If  we  combine  the  general  relation  in  (42),  for  the 
characteristic  function  of  processor  output  y,  with  the 
characteristic  function  in  (120)  for  the  sum  S  in  this  special 
case  6,  we  obtain  the  form 


fyU)  -  (1  -  i2l)K(hl{~1) 


M  K 

tt  r 

m=l  k*l( 


1  -  i2t(l  +  g  Xk  r(™>) 


-h 


exp 


il  2E 

No  m=l  k*l 


(m)' 


1  -  i2?(l  +  I!xk  r(m)) 


(151/ 


Now,  from  (67)  and  (100),  we  find  that  Kxk  covariance  matrix 

C(m)  -  r(m)  C  ,  where  C  -  [Ru ( tR-t j * fk"f  j  )  ] *  »  (152) 

and  r^^  *  1.  Then  (101)  yields 

Q(m)  -  Q  /  A(m)  -  r(m)  A  ,  (153) 

where  Q  and  A  are  the  solutions  to  the  single  Kxk  matrix  equation 

C  Q  =  Q  A  .  (154) 

Now,  define  normalized  covariance  matrix  C  as  in  (140),  with 
corresponding  matrix  equation  and  eigenvalue  matrix  A  as  in 
(141).  Then  the  eigenvalues  are  related  according  to 

Xk  *  4  Rn(°'0>  '  (155) 

leading  to  result 
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2E  * 

N  Ak 
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r (m)  =  ^  Rn(0,0)  Xk  r(m)  = 


2E 
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11  x  r(m) 
_  -k  r 


2E 


lm 


N 


^k 


(156) 


upon  use  of  (38),  (52),  and  (67).  The  quantity  E^m  is  the 
average  random  received  signal  energy  in  the  m-th  component  of 
any  one  pulse.  Also,  from  (135)  and  (137), 


2|  (m)2 

N  Ek 
o 


(157) 


since  normalized  modal  matrix  Q  is  independent  of  component 
number  m;  see  (153).  Use  of  (156)  and  (157)  in  (151)  yields  the 
characteristic  function  of  processor  output  y  in  the  desired  form 


fYU) 


=  (1  -  i20K(l5M'1) 


M  K 


m=l  k=l 


2E. 


i  -  m(i  ♦  -gis  xk) 


-h 


x  exp 


M  K 
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F(m) 
— k 


m: 


_ _  C _ _  JTf 

k=1  1  -  i2l[l  +  Xk) 


(158) 


This  result  is  a  slight  simplification  of  (136)  which  allowed 
uncorrelated  components  of  arbitrary  covariances;  thus, 
eigenvalues  (Xk)  and  eigenvectors  (Vk)  are  independent  of  m  here. 
Furthermore,  (158)  simplifies  to  the  result  in  (144)  when  the 
covariances  of  the  random  fading  components  (gm(t,f)}  in  (45)  are 
identical;  see  (64). 

The  cumulants  (Xy(P) 1  of  processor  output  y  are  found  by 
expanding  the  logarithm  of  (158)  in  a  power  series  in  i£: 
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X_^(P) 


(P-1 ) i  2P-1 


M  K  , 

EZIZ 

m=l  k=l ' 


2E 


1  + 


lm 


N 


+  p 


H  K  , 
m*=l  k=l ' 


2E 


1  + 


lm 


N 


P-1  e(rn,2 

-k 


-  K(M-2 )  for  p  1  1 


(159) 


SPECIAL  CASE  7 :  Zero  means  and  uncorrelated  components  with 
proportional  covariances 


When  we  set  means  (cm(t,f)}  to  zero,  the  characteristic 
function  in  (158)  simplifies  to 


*Y(S)  * 


(1  -  i2l) 


K(JjM-l) 


M  K 

TT  TT 

m=l  k=l 


+ 


.  (160) 


For  purposes  of  review,  the  quantity  E^  is  the  average 
random  received  signal  energy  in  the  m-th  component  of  the  k-th 
pulse,  while  is  the  received  deterministic  signal  energy  in 
the  m-th  component  of  the  k-th  pulse;  see  (50)  and  (52)  as  well 
as  fading  model  (45)  -  (47).  Eigenvalues  (Ak)  correspond  to 
normalized  covariance  matrix  C  in  (140),  while  coefficients  lv^jl 
used  in  (157)  are  the  corresponding  eigenvectors'  components;  see 
the  Kxk  matrix  equation  in  (141),  which  yields  eigenvalue  matrix 
A  and  normalized  modal  matrix  Q. 
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In  the  sequel  to  (42),  it  was  noted  that  the  more  general 
conditional  characteristic  function  of  y  in  (27)  would  be  treated 
once  the  fading  model  had  been  described.  This  latter  situation 
with  general  receiver  weights  {A^}  and  transmitted  signal 
energies  {E^}  is  considered  in  appendix  C.  In  particular,  the 
unconditional  characteristic  function  of  processor  output  y  is 
derived  for  general  correlated  fading  components  (gm(t,f))  and 
then  specialized  to  the  case  of  uncorrelated  fading  components. 
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EXAMPLES  OF  PROBABILITY  DENSITY  OF  POWER-SCALING  q(t,f) 

In  an  earlier  section,  the  characteristic  function  of 
power-scaling  variate  q  =  q(t,f)  was  derived  for  fading  model 
(45),  yielding  general  result  (76).  This  result  was  then 
simplified  for  the  seven  special  cases  noted  there.  In  this 
section,  we  will  present  some  numerical  examples  of  typical 
probability  density  functions  of  q  and  thereby  illustrate  the 
generality  and  variety  of  fading  model  (45). 

The  only  case  we  will  consider  in  this  section  is  where  the  M 
random  components  {gm(t,f)J  in  (45)  are  uncorrelated  with  each 
other  and  have  proportional  covariances.  This  is  special  case  6; 
see  (67).  Also,  for  notational  convenience,  we  consider  the 
normalized  (power)  random  variable 


♦  = 


q _ 

Ru(0,0) 


(161) 


where 


3  =  **(0,0)  =  Rii(0,0)  r(m)  for  1  i  m  S  M  .  (162) 

m  3m '  mm  ''ll'' 


Here,  we  used  (47)  and  (67).  Then,  from  (86),  the  characteristic 
function  of  random  variable  $  is 


V51  *  fq 


exp 


M 


H 


2  (m) 

1_  r '  ' 
'm 


m=T  l-i2Sr(m) 

(163) 


where  we  defined  (dimensionless)  normalization  constants 
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c  C_ ( t , f ) 

.  m  _m  m  nr  ' 

"  °m  ’W°'0>‘ 


for  1  *  m  i  M  . 


(164) 


For  M  *  1,  the  general  result  in  (163)  simplifies  to 


f.U)  -  (1  -  i20"*  exp 


u  n2 

1  -  i25 


(165) 


for  which  the  corresponding  PDF  (probability  density  function)  is 
p^(u)  -  ( 2nu ) ”**  exp|-  ^[u  +  t)2)]  coshfl^/u)  for  u  >  0  .  (166) 

J- 

The  amplitude-scaling  variate,  r^  »  q£  from  (43),  has  the 
corresponding  PDF  for  normalized  version 


0  “  "  [R^foToll  " 


(167) 


namely 


Po ( u )  -  2  u  p.(u  ) 


exp|-  ^[u2  +  ))2|J  coshthjU)  for  u  >  0  . 


(168) 


This  PDF  is  displayed  in  figure  2  for  ^  -  0 ( -5)3.  At  the  origin 

u  2 

this  PDF  is  finite  and  nonzero,  with  value  ( 2/n )  exp(~h^/2), 

indicating  the  possibility  of  occasional  deep  fades  (unless  is 
large).  As  parameter  ^  gets  large,  this  PDF  approaches 
Gaussian;  in  fact,  directly  from  (168)  and  figure  2,  we  obtain 


p0(u)  ~  (2JI)"*5  exp[-  ^(u  -  hi)2]  for  \\1  >  3  .  (169) 
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The  corresponding  PDF  for  normalized  power-scaling  variate 

2 

$  =  q/Rn(0,0)  =  q /aj  was  given  in  (166)  and  is  plotted  in  figure 

3.  It  is  seen  to  possess  a  substantial  cusp  at  the  origin,  which 

l  2 

behaves  as  (2nu)  exp( -1)^/2).  This  is  the  most  severe  case  of 
deep  fading  that  model  (45)  can  yield,  namely  for  M  =  1 . 

When  M  =  2,  that  is,  two  components  in  fading  model  (45),  the 
characteristic  function  in  (163)  reduces  to 

v<>  - 11  - i2£)'”  l1  - i2£  ^  -4^4 ♦  • 

(170) 

However,  instead  of  Fourier  transforming  this  result,  the  PDF 

of  $  is  best  found  by  reverting  to  definition  (45),  namely 
2  2 

q  ■  (Cj  +  gj)  +  (c2  +  g2)  ,  and  performing  the  probability 
integrals  directly  in  the  9^*92  plane-  The  result  is 

n 

P*(u)  *  ^  j  ?2)j^  J  dt  exP[‘  l[c2  +  s2) ]  for  u  >  0  ,  (171) 

where  auxiliary  functions 

C  s  u1*  cos(t)  -  /  S  s  ^  sin(t)  -  h2  •  (172) 

As  a  special  subcase  here,  for  M  =  2,  let  r^  =  1;  that  is, 
let  both  random  components  (gm(t,f)}  have  equal  power.  Then 
(170)  yields  characteristic  function 
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2 

%<*>  *  (1  -  i2trl  exP[l  -  ktl  ' 

where  constant 

»2  -  *i  +  "j  -  siTTs^y  £°r  m  -  2  . 

The  corresponding  probability  density  function  for  4  is 
p^(u)  -  |  IQ(h^u)  exp[-  ^(u  +  h2) ]  for  u  >  0  . 

That  for  0  «  follows  immediately  as 

pQ(u)  =  u  IQ(hu)  exp[-  ^(u2  +  h2)]  for  u  >  0  . 

The  PDF  in  (175)  is  displayed  in  figure  4.  Now,  the  origin 


(173) 

(174) 

(175) 

(176) 

value 


TR  10041 


2 

is  finite,  namely  .5  exp(-^  /2 ) ,  in  contrast  with  figure  3  for 
M  =  1 .  The  cusps  that  were  in  figure  3  are  now  absent  for  M  =  2. 

When  M  =  2  but  r^2^  ^  1,  a  closed  form  like  (175)  is  not 

available  for  the  PDF  of  Instead,  we  must  revert  to  the 

general  relation  in  (171)  -  (172)  and  perform  the  numerical 
integration  for  each  set  of  parameter  values  of  interest.  An 
example  for  r*2^  *  1/2  and  *)2  =  ^1  is  9iven  in  figure  5,  and  a 

second  example  for  r*2*  =  1/2  and  h2  *  h 2  is  given  in  figure  6. 

A  common  feature  of  all  these  densities  is  the  finite  nonzero 
values  at  the  origin.  That  is,  deep  fades  are  still  common  for 
fading  model  (45)  when  M  =  2. 

For  larger  M  and  completely  arbitrary  parameters,  it  is 
necessary  to  numerically  Fourier  transform  general  characteristic 
function  (163)  in  order  to  determine  the  PDF  of  +.  However,  for 
the  special  case  where 

M  is  even  and  r(m)  #  1  for  1  £  m  £  M  ,  (177) 

then  (163)  reduces  to 

V51  *  (1  -  i2«)'M/2  ‘HZ-  ml  < 

where  constant 

0  M  9 

H2  -  El  "l  ■ 

m=l 

The  corresponding  PDF  of  $  is  then 

p+(u)  =  \  (^)N  IN(h^u)  exp [-  -|(u  +  h2)  ]  for  u  >  0  ,  (180) 


(178) 


(179) 
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Figure  6.  PDF  of  q /o^  for  M  =  2,  r^^  *  1/2,  h 2  *  hj_/2 
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where  N  =*  M/2  -  1.  This  reduces  to  (175)  for  M  *  2. 

The  PDF  for  0  =  /$  follows  as  given  by  the  rule  given  in 
(168),  namely 

N 

Pe(u)  «  u  (^)  IN(hu)  exp  j[u2  +  h2)]  for  u  >  0  .  (181) 

This  reduces  to  (176)  for  M  =  2. 

Two  examples  of  PDF  (180)  for  random  variable  $  are  displayed 
in  figure  7  for  M  *  4  and  in  figure  8  for  M  =  6,  respectively. 

The  former  PDF  is  zero  at  the  origin,  while  the  latter  is  zero 
and  has  zero  first  derivative  at  the  origin.  Thus,  deep  fades 
are  less  likely  for  the  larger  values  of  M  in  fading  model  (45). 

A  general  procedure  for  determining  the  probability  density 
function  p^(u)  of  random  variable  $  from  characteristic  function 
f^(£)  in  (163),  for  arbitrary  M,  { r ^ m ^ is  presented  at 
the  end  of  appendix  D.  It  employs  some  efficient  recursions  for 
fast  and  accurate  numerical  evaluation  of  p^(u). 
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PROBABILITY  DISTRIBUTION  OF  PROCESSOR  OUTPUT  y 

The  characteristic  functions  fy(S)  of  processor  output  y, 
under  a  variety  of  special  cases,  were  derived  in  an  earlier 
section;  in  particular,  see  (124)  -  (160).  Although  these  are 
compact  closed  forms,  they  can  encounter  computational  problems 
if  we  attempt  to  directly  utilize  Fourier  transforms  to  determine 
the  corresponding  probability  density  functions  and  exceedance 
distribution  functions.  In  particular,  since  the  characteristic 
functions  decay  only  as  (  as  (  ■»  ®,  truncation  error  can  become 
a  significant  problem,  especially  for  small  K,  the  number  of 
signal  pulses. 


EXPANSION  FOR  DISTRIBUTION  FUNCTION 

In  appendix  D,  a  series  expansion  for  and  the 

corresponding  probability  density  function  p^(u)  are  derived  in  a 
form  which  involves  only  positive  expansion  coefficients  and 
terms.  Furthermore,  efficient  and  accurate  recursions  are 
developed  for  the  evaluation  of  the  coefficients,  the  probability 
density  function,  and  the  exceedance  distribution  function. 

As  an  example  of  the  general  procedure,  the  results  for 
general  characteristic  functions  (136)  and  (158)  are  summarized 
here.  Define,  for  1  £  m  S  M,  1  £  k  £  K, 


mk 


^Elm  i (m) 

nr  or 


2E 


lm 


mk 


N 


-k  ' 


e  =  E{ra) 
mk  -k 


(182) 


63 


TR  10041 


Then,  the  exceedance  distribution  function  of  processor  output  y 
is  given  by  (D-16)  as 


w 

Prob ( y  >  u>  -  P  ZZ  gp  for  u  >  0  , 


(183) 


where  the  quantities  required  are  found  in  the  following  fashion. 

■h 

,  (184) 


H 

K,  \ 

F  = 

n  t 

m=l  k 

I[(l  +  Tmk) 

a 


o 


'mk 


1  +  T 


mk 


i 


a  = 


1  >M  )K 


PP-1 

mk 


m-l  k=l  (1  +  Tmk) 


'mk 


1  +  T 


mk 


for  p  1  1 


( 185 


90  -  exp(aD)  ,  9p  4  C  "  °n  Vn  £or  P  *  1  '  <186> 

n=  1 

n  k 

H  (x)  =  exp(-x)  )  ~  £7  for  n  1  0  ,  x  1  0  .  (187) 

n  k=0  K* 

Since  the  eigenvalues  of  a  covariance  matrix  can  never  be 
negative,  it  is  seen  from  (185)  that  all  the  (Op)  are 
nonnegative,  with  the  exception  of  aQ.  However,  aQ  is  used  only 
once  in  (186)  to  generate  a  positive  gQ,  while  the  recursion  in 
(186)  for  (gp)  utilizes  only  nonnegative  quantities.  Thus,  all 
the  coefficients  (gp)  and  scale  factor  F  in  (183)  are 
nonnegative.  The  function  Hn(x)  in  (187)  is  obviously  positive 
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for  x  1  0,  thereby  guaranteeing  that  no  negative  terms  will 
appear  in  expansion  (183).  Furthermore,  there  is  a  very 
efficient  recursion  for  the  (Hn(x)};  see  (D-18).  A  sample 
program  for  the  evaluation  of  a  modified  form  of  (183)  is  given 
in  appendix  D. 

In  the  absence  of  signal,  all  the  {TmJc}  and  {emkl  in  (182) 
are  zero.  Then,  all  the  {<Xp}  and  ( g^ }  become  zero  except  for  gQ 
which  is  1.  Expansion  (183)  for  the  exceedance  distribution 
function  of  processor  output  y  then  reduces  to  simply  the  one 
term  HRl(u/2),  which  is  consistent  with  the  noise-only 
characteristic  function  (1  -  i2£)  for  processor  output  y;  see 
(42),  for  example.  That  is,  the  false  alarm  probability  is 


An  error  bound  for  sum  (183),  terminated  at  the  p  =  N  term, 
is  given  in  (D-20).  Also,  the  modifications  required  to  treat 
the  slightly  different  forms  of  characteristic  functions 
encountered  in  (124)  and  (144)  are  presented  in  (D-23)  -  (D-26). 
Finally,  a  refinement  of  the  expansion  procedure  in  (183),  which 
is  more  rapidly  convergent  and  useful  for  larger  signal-to-noise 
ratios,  is  given  in  (D-27)  -  (D-36). 
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FUNDAMENTAL  INPUTS  AND  COMPUTATIONAL  PROCEDURE 

In  the  remainder  of  this  section,  we  will  restrict  attention 
to  characteristic  function  (158)  which  pertains  to  special  case 
6,  namely  uncorrelated  components  with  proportional  covariances; 
see  (67).  However,  before  we  list  the  fundamental  inputs  that 
are  required  to  conduct  the  numerical  evaluation  of  the 
exceedance  distribution  function  of  processor  output  y  in  this 
case,  we  make  an  additional  modification  for  later  convenience  in 
plotting  and  comparison. 

From  (52),  (38),  and  (67),  we  have 

Eki»  ’  Elm  -  E  Rll<°'°>  r<m)  •  US’) 

Then,  using  (53),  the  average  random  received  signal  energy  in 
one  pulse  is 

M  M 

E1  -  EZ  E,m  -  E  Ru(0,0)  HZ  r'm'  .  (190) 

m= 1  m=l 

This  enables  us  to  express  the  signal-to-noise  ratio  measure  as 
Elm  El  (m) 

T  ‘  r  +.  ■  where  »„  11  H  —  for  1  S  m  S  M  .  (191) 

°  III  r(m) 

m=l 

The  quantities  represent  the  fractional  strengths  of  each  of 

the  M  random  components  (g  (t,f)}  in  fading  model  (45). 

The  fundamental  inputs  required  to  evaluate  the  exceedance 
distribution  function  of  processor  output  y  are  now 
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K,  number  of  signal  pulses  in  figure  1, 

M,  number  of  fading  components  in  (45) , 

E1/Nq/  signal  energy  to  noise  density  ratio, 
lr^m) }  for  1  £  m  £  M;  see  (67), 
lDkm/Nol  for  1  *  k  *  K/  ISmSM;  see  (50), 
Ru(t,v)/R11(0,0);  see  (47), 

i tR J  and  { f k }  for  1  i  k  1  K;  see  (44).  (192) 

The  quantity  E1/Nq  is  a  measure  of  the  receiver  input  average 
random  signal-to-noise  ratio  for  one  signal  pulse,  while  r^ra^  is 
the  relative  strength  of  the  m-th  fading  component.  Ratio 
is  a  measure  of  the  receiver  input  deterministic  signal-to-noise 
ratio  for  the  k-th  signal  pulse  and  m-th  component.  The  function 
Rll(T 'v)/Rii(0,0)  is  the  normalized  fading  covariance  function 
for  the  medium  at  time  separation  t  and  frequency  separation  v. 
Parameters  t^  and  f^  are  the  time  and  frequency  locations, 
respectively,  of  the  k-th  signal  pulse  in  time, frequency  space. 

The  first  quantity  that  must  be  computed  is  the  normalized 
Kxk  covariance  matrix  C  given  by  (140).  Then,  its  eigenvalue 
matrix  A  and  normalized  modal  matrix  Q  are  found  by  solving 
characteristic-value  matrix  equation  (141).  This  yields 
eigenvalues  and  eigenvectors  {V^j  for  1  £  k  £  K.  The 

components  of  vector  are  v^,  •  •  •  ,v^K;  see  (102).  We  then 
compute  from  (142)  for  1  £  k  £  K,  1  £  m  £  M.  At  this  point, 

the  parameters  {Tmk)  and  {e^j  in  (182)  can  be  computed  for 
1  £  k  £  K,  1  £  m  £  M,  and  the  procedure  in  (183)  -  (187)  can  be 
employed.  In  particular,  (182)  and  (191)  yield 
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2E. 

Tmk  "  IT  ’♦m  -k  for  1  *  m  *  M  #  1  £  k  £  K  .  (193) 

o 

NUMERICAL  PERFORMANCE  RESULTS 

The  first  set  of  examples  are  selected  to  enable  a  comparison 

with  the  approximate  procedure  and  results  in  [12].  In  figure  9, 

below,  we  plot  the  required  per-pulse  input  SNR  (signal-to-noise 

ratio)  measure  E1/Nq  (in  dB)  versus  the  number  of  signal  pulses 

K,  for  values  of  the  adjacent-pulse  normalized  (amplitude) 

covariance,  Cov,  equal  to  0,  .5,  /.5,  and  1.  In  particular,  the 

covariance  R^^(t,v)  was  taken  as  exponential  in  t,  and  the  signal 

pulses  have  no  frequency  shifts  (f^J  and  are  equally  spaced  in 

time  locations  [tkj.  Also,  the  KM  deterministic  signal-to-noise 

ratios  [Dj^/N  }  are  zero,  and  the  M  strengths  {r^  }  are  all 

equal  to  1;  this  duplicates  the  situation  in  [12].  The 

particular  case  in  figure  9  pertains  to  false  alarm  probability 

Pp  »  IE-6 ,  detection  probability  PD  =  .5,  and  M  =  1  fading 

component;  this  last  choice  corresponds  to  parameter  m  =  1/2  in 

[12].  The  results  here  in  figure  9  for  Cov  =  0  and  Cov  =  1  agree 

precisely  with  p  =  0  and  p  =  1  in  [12;  figure  11],  as  expected, 

while  the  results  for  Cov  *  /.5  are  in  rather  good  agreement  with 

those  for  p  ®  .5.  This  selection  of  parameters  reflects  the 

2  . 

property  that  p  =  Cov  ;  see  the  section  entitled  Program  m 
appendix  D  for  additional  details. 

The  only  change  in  figure  10  is  to  require  a  larger  detection 
probability,  namely  PD  *  .9.  Although  the  Cov  *  0  and  1  results 
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Figure  9.  Required  SNR  for  Pp  «  IE-6,  PD  *  .5,  M  =  1 


69 


TR  10041 


agree  with  [12;  figure  12],  the  exact  results  here  for  Cov  «  /.5 
are  markedly  different  from  the  approximate  results  for  p  =  .5  in 

the  earlier  work.  For  example,  2.6  dB  less  is  required  for  K  =  2 

and  K  a  3,  while  2.3  dB  less  is  required  for  K  »  4.  That  is,  the 
earlier  approximation  in  [12]  was  somewhat  pessimistic  in  its 
performance  predictions. 

We  must  observe  from  figure  10  that  the  effect  of  correlated 
fading  is  not  overly  significant  until  the  normalized  covariance 
approaches  1.  For  example,  for  K  =  2,  the  cost  of  Cov  increasing 
from  0  to  /.5  is  1.4  dB,  while  the  cost  of  Cov  going  from  /.5  to 
1  is  an  additional  6.7  dB.  The  distinction  is  even  greater  for 

K  »  10,  with  the  same  comparison  requiring  1.6  dB  versus  an 

additional  12.9  dB. 

The  example  in  figure  11  pertains  to  M  *  2  fading  components, 
(which  corresponds  to  m  =  1  in  [12]);  all  other  parameters  are 
the  same  as  figure  9  above.  These  results  can  be  compared 
directly  with  [12;  figure  7];  they  reveal  identical  performance 
for  Cov  88  0  and  1,  and  rather  good  agreement  for  Cov  «  /.5 
(versus  p  *  . 5 ) . 

When  the  detection  probability  PD  is  increased  to  .9  in 
figure  12,  these  exact  results  reveal  that  lower  values  of  E1/Nq 
are  required  than  the  approximation  in  [12;  figure  8]  predicted, 
at  least  for  Cov  =  /.5  (p  =  .5).  Also,  as  was  seen  in  figure  10, 
the  cost  of  normalized  covariance  Cov  approaching  1  is  very 
significant  in  terms  of  increased  signal  level;  that  is,  the 
Cov  =  /.5  curves  are  well  below  the  Cov  =  1  curves  in  figures  10 
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Figure  11.  Required  SNR  for  Pp  =  IE-6,  PD  -  .5,  M  *  2 
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and  12,  especially  for  the  larger  K  values. 

The  detrimental  effect  of  highly  correlated  fading  pulses  is 
studied  quantitatively  in  figure  13,  where  normalized  covariance 
Cov  is  varied  from  0  to  1  as  the  other  various  parameters  are 
kept  fixed.  The  number  K  of  signal  pulses  is  kept  at  2.  These 
plots  reveal  that  for  detection  probability  PD  =  .5,  increased 
Cov  does  not  lead  to  a  significantly  large  increase  in  required 
input  signal  level.  However,  for  the  better  detection 
probability  of  PD  =  .9,  higher  covariances  can  be  very  damaging, 
requiring  additional  signal  strength  to  maintain  the  desired 
level  of  performance.  For  example,  as  Cov  increases  from  .9  to  1 
in  the  uppermost  example  in  figure  13,  the  signal  must  be 
increased  by  4.9  dB.  (The  kinks  in  the  curves  are  due  to 
discretization  of  the  abscissa  at  increment  .1  for  Cov.) 

When  K  is  increased  to  6,  the  results  in  figure  14  reveal 
this  effect  in  a  more  pronounced  fashion.  An  additional  9  dB  is 
now  required  when  Cov  is  increased  from  .9  to  1  for  the  upper 
curve. 

All  of  the  above  results  have  had  deterministic  signal-to- 
noise  ratio  measures  (Dkm/NoJ  equal  to  zero.  The  program  listed 
in  appendix  D  has  the  capability  of  incorporating  arbitrary 
values  for  these  parameters  as  well  as  others,  such  as  { r ^ ^ } . 

An  example  where  all  the  potential  is  exercised,  and  all 
parameters  have  nonzero  values,  is  displayed  in  figure  15.  Here, 
the  detection  probability  PD  is  varied  from  .5  to  .999  and  the 
required  E^N  is  calculated  (in  dB).  It  is  seen  that  a  sharp 


72 


o  ,1  ,1  ,3  .+  _  .S'  .<  J  ,»  .1  I 

Lev 

Figure  14.  Required  SNR  versus  Cov,  K  =  6 
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increase  is  observed  as  PD  increases  above  .9,  eventually  tending 
to  ®  as  PD  -*  1.  The  particular  parameter  values  are  listed 
below: 


PF  =  IE-6,  K  =  4,  M  =  2, 

r(m)  =  i  for  1  $  m  i  M, 

Dkm/No  *  1  for  1  *  k  S  K,  1  1  m  S  M, 

|tk|  =  1,2, 3,4,  {fk}  =  1,4, 2, 3, 

2  2 

Cov  =  exp(-  XT  -  Yo)  - 


(194) 


Figure  15.  Required  SNR  versus  PD 
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SUMMARY 

The  characteristic  function  of  processor  output  y  in  figure  1 
has  been  derived  in  closed  form  for  a  wide  variety  of  fading 
conditions  and  signal  formats.  The  corresponding  probability 
density  functions  and  exceedance  distribution  functions  of  y  have 
then  been  expanded  in  convergent  series,  by  means  of  a  novel 
expansion  technique  which  make  efficient  use  of  recursions  for 
rapid  and  accurate  evaluation.  These  series  typically  require  on 
the  order  of  30  to  60  terms  for  ten  decimal  accuracy.  A  program 
is  given  which  incorporates  all  these  features,  including 
allowance  for  both  deterministic  and  random  components  of 
arbitrary  strengths  in  the  fading  medium. 

One  of  the  most  useful  results  pointed  out  by  this  study  is 
that  the  degree  of  correlation  between  the  fading  signal  pulses 
can  be  fairly  significant  without  suffering  great  degradations. 
That  is,  when  the  normalized  covariance  approaches  1,  meaning 
that  all  signal  pulses  tend  to  fade  together,  the  performance 
losses  are  potentially  large;  however,  for  coefficients  below  .5, 
the  losses  are  not  too  significant. 

Comparison  of  these  exact  results  with  an  earlier  approximate 
procedure  [12]  reveals  that  the  earlier  approach  generally  gives 
pessimistic  predictions  of  performance  when  the  normalized 
covariance  of  the  fading  is  intermediate  between  0  and  1.  In 
some  cases,  the  discrepancy  can  be  several  dB.  This  result 
indicates  and  emphasizes  the  need  for  accurate  treatment  of 
systems  which  must  perform  well,  that  is,  yield  high  detection 
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probabilities  in  fading  media.  An  extension  of  this  work  to  a 
fading  medium  in  which  the  background  noise  level  is  unknown  and 
must  be  estimated  from  a  finite  sample  in  a  noise-only  region  of 
time, frequency  space,  is  currently  underway  by  this  author  and 
will  be  reported  on  shortly. 

The  major  result  utilized  here  is  the  characteristic  function 
of  a  quadratic  form  in  correlated  nonzero-mean  Gaussian  random 
variables;  this  result  is  presented  in  appendix  B.  It  relies  on 
the  ability  to  solve  the  generalized  eigenvalue  problem;  this 
latter  procedure  and  solution  is  presented  in  appendix  E.  A 
related  problem  involving  a  slightly  more  general  bilinear  form 
is  treated  in  appendix  F.  Finally,  the  characteristic  function 
of  the  most  general  complex  form  with  both  first-order  and 
second-order  terras  is  solved  in  appendix  G.  These  results  are 
not  utilized  in  this  technical  report  but  are  presented  for 
completeness  and  for  possible  future  use  and  reference. 
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APPENDIX  A.  ANALYTIC  AND  COMPLEX  ENVELOPE  PROCESSES 

Let  n(t)  be  a  stationary  zero-mean  real  random  process  with 
covariance  Rn(r)  and  double-sided  spectrum  Gn(f): 

Rn(t)  =  n(t)  n(t-r)  ,  Gn(f)  =  J  dT  exp(-i2nfr)  Rn(t)  .  (A-l) 

The  analytic  process  n+(t)  is  generated  by  eliminating  the 
negative  frequencies  in  n(t)  and  by  doubling  its  positive 
frequency  contributions;  that  is,  n(t)  is  passed  through  a  filter 
with  transfer  function  2  U(f),  where  U  is  the  unit  step  function. 


Figure  A-l.  Generation  of  Analytic  Process 

The  analytic  process  can  be  expressed  in  terms  of  its  real 
and  imaginary  parts  according  to 

n+( t )  =  n(t)  +  i  nH(t)  ,  (A-2) 

where  nH(t)  is  the  Hilbert  transform  of  n(t).  The  complex 
envelope  process  n(t)  is  obtained  by  frequency  down-shifting 
analytic  process  n+(t),  in  order  to  center  its  one-sided  spectrum 
about  f  =  0.  Thus,  if  fQ  is  a  representative  center  frequency  of 
n+ ( t ) ,  then 

n(t)  =  n+(t)  exp(-i2nfQt)  .  (A-3) 
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The  original  real  process  n(t)  then  follows  from  (A-2)  and  (A-3): 

n(t)  =  Re|n+(t)j  =  Re|n(t)  exp(i2nfQt)j  .  (A-4) 

The  spectrum  of  n+(t)  follows  from  figure  A-l  as 

Gn+(f>  =  4  u(f)  Gn(f)  '  (A-5) 

while  the  spectrum  of  n(t)  is  obtained  from  (A-3)  according  to 

G  ( f )  =  G  (f+f  )  =  4  U(f+f  )  G  ( f+f )  •  ( A-6 ) 

n  n+  o  '  o  n'  o 

The  covariances  corresponding  to  (A-5)  and  (A-6)  are 

n+(t)  n* ( t-r )  and  n(t)  n*(t-r)  ,  (A-7) 

respectively.  The  two  complementary  covariances  are  zero;  that 
is, 

n+(t)  n+(t-x)  =  0  ,  n(t)  n(t-t)  =  0  .  (A-8) 

These  results  follow  from  the  fact  that  transfer  function  2  U(f) 
in  figure  A-l  is  single-sided;  that  is,  U(f)  =  0  for  f  <  0. 

A  situation  that  frequently  arises  in  practice  is  where  n(t) 
is  a  noise  process  with  a  spectrum  Gn(f)  that  is  flat  in  a  broad 
band  of  width  W  in  the  neighborhood  of  fQ,  which  essentially 
covers  the  signals  and  filters  of  interest.  The  spectra  of  the 
various  processes  are  illustrated  in  figure  A-2,  where  is  the 
double-sided  noise  spectral  density  level  of  n(t)  in  the 
neighborhood  of  The  spectrum  Gn(f)  of  the  complex  envelope 

n(t)  is  flat  in  the  neighborhood  of  f  =  0  and  has  level  4 
watts/Hz.  Therefore,  its  covariance 
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Figure  A-2.  Spectra  of  Analytic  and  Complex  Envelope  Processes 
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Rn*T)  =  J  df  exP(i2nft)  Gn(f)  (A-9) 

is  (relative  to  the  signal  and  filter  time  functions)  a  sharp 
pulse  centered  at  r  -  0,  with  area  Gn(0)  =  4  N^.  Then,  a  good 
approximation  for  many  purposes  is  to  say  that 

Rn(t)  •  4  Nd  S(t)  •  (A- 10 ) 

There  is  no  need  to  let  W  and  f  tend  to  infinity  in  order  to 
utilize  this  result;  it  only  requires  that  W  be  somewhat  larger 
than  the  bandwidths  of  the  signals  and  filters  of  interest. 

In  terms  of  the  one-sided  noise  spectral  density  level  of 
n(t),  that  is,  NQ  *  2  Nd,  approximation  (A-10)  becomes 

Rn(T)  5  2  No  &(T)  *  (A_11) 

This  result  may  be  compared  with  [19;  page  48,  (3.10)  -  (3.11) 
and  page  72,  (6.21)  -  (6.22)],  where  his  N  is  our  Nq. 

It  should  be  noted  that  original  covariance  Rn(r)  in  (A-l) 
cannot  be  recovered  from  this  approximation.  That  is,  via  (A-4), 

Rn(0  *iRe(Rn(T)  exP(i2nfor))  =  (A-12) 

■  j  Re(4  Nd  6 ( r )  exp(i2nfor)j  =  2  Nd  6(t)  ,  (A-13) 

which  is  incorrect.  The  flaw  is  that  the  use  of  6(t),  which  is 

tantamount  to  W  ->  ®,  must  be  accompanied  by  having  let  fQ  -»  ®,  in 

addition.  That  is,  exp(i2nfQr)  in  (A-12)  must  vary  as  fast  as 

R  (t)  in  order  for  the  narrowband  representation  to  be  valid, 
n ' 

The  correct  end  result  for  fQ  =  W/2  is  Rn(T)  &(T)  as  W  -*  ®. 
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APPENDIX  B.  CHARACTERISTIC  FUNCTION  OF  QUADRATIC  AND  LINEAR 
FORM  IN  CORRELATED  NONZERO-MEAN  GAUSSIAN  RANDOM  VARIABLES 

The  Nxl  real  vector  X  is  composed  of  correlated  Gaussian 
components  with  mean  vector  E  and  covariance  matrix  C;  that  is, 

X  =  E  ,  Cov(X)  ={X-E)(X-E)T=C.  (B-l) 

Covariance  matrix  C  is  Nxn,  real,  symmetric,  and  nonnegative 
definite.  Real  Nxl  vector  E  is  completely  arbitrary. 

We  shall  be  interested  in  obtaining  the  characteristic 
function  of  the  quadratic  and  linear  real  form  in  X  given  by 

q  =  q( X )  =  XT  B  X  +  2  AT  X  ,  (B-2) 

where  nxn  matrix  B  is  real,  symmetric,  and  positive  definite, 
while  real  Nxl  vector  A  is  completely  arbitrary.  By  completing 
the  square,  (B-2)  can  be  written  in  the  alternative  form 

q  =  (X  +  B-1  A)T  B  (X  +  B~ 1  A)  -  AT  B-1  A  .  (B-3) 

T  -1 

Therefore,  the  minimum  possible  value  of  q  is  -  A  B  A.  The 
case  of  an  indefinite  matrix  B  is  undertaken  below  (B-42). 

Random  vector  X  is  composed  of  correlated  Gaussian  random 
variables;  its  probability  density  function  is  [20;  section  8-3] 

p(X)  =  (2n)'N/2  (det  C)"1/2  exp[-  ^(X  -  E)T  C_1  (X  -  E)j.  (B-4) 

The  characteristic  function  of  interest  is  then  given  by 

fq[l)  *  exp(Hq)  -  J  dX  p(X)  exp[Hq(X)]  .  (B-5) 
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For  convenience,  we  introduce  variable 

z  =  i2l  .  ( B-6 ) 

Substitution  of  (B-2)  and  (B-4)  in  (B-5)  yields 


fqU)  =  (2n)'N/2  (det  C)"1/2  x 
x  J  dX  exp[-  ^{X  -  E)T  C'1  (X  -  E)  +  jz  XT  B  X  +  z  AT  x]  = 


*  (2n)"N/2  (det  C)"1/2  x 


X  J  dX  exp|-  ixT(C_1  -  zB)X  +  (C-1  E  +  zA)TX  -  ^ET  C_1  e]  .  (B-7) 


Now  we  have  the  result  (20;  section  8-3] 

J  dX  exp|-  |xT  u  X  +  VT  x]  =  (2n)N/2  (det  U)“1/2  exp[^VT  U_1  v]  . 

(B-8) 

This  enables  the  reduction  of  the  integral  in  (B-7)  to 

fq(f;)  -  (det  C)-1/2  (det(C-1  -  zB))"1/2  x 

[1-1  T  - 1  -1-1  IT  -11 

i(C  E  +  zA)  (C  -  zB)  (C  E  +  ZA)  -  ^E  C  El  .  (B-9) 

Although  closed  form,  (B-9)  is  not  too  useful  numerically  because 
it  requires  an  inverse  of  matrix  C-*  -  zB  for  each  new  value  of  z 
(=  i2£)  of  interest. 

A  much  more  compact  and  useful  form  of  (B-9)  can  be  obtained 
by  means  of  the  following  procedure.  For  the  given  covariance 
matrix  C  and  quadratic- form  matrix  B,  solve  the  generalized 
characteristic-value  equation 
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C  Q  -  B  1  Q  A  ( B-10 ) 

for  Nxn  eigenvalue  matrix  A  and  normalized  modal  matrix  Q,  where 

A  =  diag[A1  ...  XN]  f  Q  -  [Vx  ...  VN]  ,  Vn  *  [vnl  ...  vnN]T  , 

(B-ll) 

and  Nx 1  column  vector  Vn  is  the  n-th  eigenvector  with  scalar 
components  { vnp I  -  A  procedure  for  obtaining  the  solutions  A  and 
Q  to  (B-10)  is  given  in  appendix  E,  when  matrix  B  is  positive 
definite. 

(It  should  be  noted  that  (B-10)  does  not  have  exactly  the 
same  solutions  as  the  equation  B  C  Q'  =  Q'  A';  in  fact,  B  C  is 
not  generally  symmetric,  even  if  B  is  diagonal  [18;  page  79, 
(253)].  The  connection  is  A'  *  A,  Q'  “  Q  D,  with  D  diagonal.) 

Then,  we  have  the  two  very  important  properties  [18;  pages 
74  -  77]  of  solutions  Q  and  A: 

QT  B_1  <2  =  1,  QTCQ  =  A.  ( B-12 ) 

By  means  of  these  two  relations,  a  number  of  simplifications  of 
( B—  9 )  are  possible.  We  begin  by  observing  that 

-1  -T-l  T  -T  -1  -1  -IT 

B  =  Q  Q  ,  B  =  Q  Q  ,  C  =  Q  A  Q  ,  C  -  Q  A  Q  .  (B-13) 

(Notice  that  QT  Q  ^  I . )  We  now  employ  (B-13)  in  (B-9)  to  get 


-1  -1  T  -1  -1  -T  -1  -l  -l 

C  -  zB  =  Q  (A  -  zl)  Q  ,  (C  -  zB)  -  Q  A  (A  -  zl)  Q  . 

(B-14 ) 


At  the  same  time,  using  (B-13), 
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det(C)  det(C-1  -  z  B)  =  det(I  -  z  B  C)  =  det(I  -  z  Q  A  Q_1 )  = 


=  det ( I  -  z 


A)  ■  Ht1  -  *  *») 


( B— 1 5 ) 


Meanwhile,  if  we  denote  twice  the  argument  of  the  exponential 
in  (B-9)  by  t,  we  have,  by  means  of  (B-14),  the  alternative 
expression 

-1  T  -T  -1  -1  -1  -1  T  -1 

t  «  (C  E  +  zA)  Q  1  (A  -  zl)  Q  (C  E  +  zA)  -  E  C  E  . 

(B-16 ) 

Now,  by  means  of  ( B— 13),  develop  the  term 

-1  -1  -1  -1  -1  -IT  T  -1 

Q  (CE+zA)=QCE+zQA=AQB+  zQ  B  A  = 


=  A  1  E  +  zA  , 


( B-17 ) 


where  we  have  defined  N*1  vectors 


E  =  QT  E  =  [tl  . 


T  T 

*  EN>  '  En  =  Vn  E  ' 


A  *  QT  B'1  a  -  [ax  ...  a^)T  ,  an  =  V*  B'1  a  .  (B-18) 


Then  t  in  (B-16)  becomes,  again  using  (B-13) 


t  =  (A*1  E  +  zA ) T  (A-1  -  zI)"A  (A-1  E  +  zA)  -  ET  Q  A  1  QT  E  = 

-1  t  - 1  -1-1  T-l 

=  (A  E  +  zA)  (A  -  Zl)  (A  E  +  zA)  -  E  A  E  = 


N  f e  /A  +  z  a  )2  N  z2 
— •  Inn  nj 

—  lA  -  z 
n=l  n 


n 


n=l  ^n 


N  e_  +  a)  N  9 

EZ  1-  z  i  -  2  S  “n  • 
n=l  n  n=l 


(B-19) 
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By  collecting  all  these  results  together,  we  can  express 
the  characteristic  function  in  (B-9)  in  the  compact  form 

*1  . ( B-20 ) 

The  last  constant  in  (B-20)  has  the  alternative  representation 

N 

Ll  9  T  T  —1  T  —1  T  —1 

-  ZZZ  a^  =  -AiA=-AiBiQQiBiA=-AABAA/  (B-21) 
n=l 

where  we  used  (B-18)  and  (B-13);  this  is  just  the  residual 
constant  encountered  in  (B-3).  Thus,  the  characteristic  function 
of  the  leading  (nonnegative)  quadratic  form  in  (B-3)  is  just 
(B-20)  without  the  second  summation  inside  the  exponential. 
Additional  relations  available  from  (B-18)  and  (B-13)  are 

n=l 

N 

”  T  T  T  -1  T  T 

}  £a=  E1  A  =  E1  Q  Q1  BA  =  EA  =  AE  .  (B-22) 

n=l 

However,  we  need  the  constants  { cxn }  and  { en }  individually  for 

characteristic  function  (B-20).  The  quantities  { 2 }  and 
2 

| ( c  +  an)  )  should  be  computed  once  and  stored  m  arrays  prior 
to  the  computation  of  fg(£)  i-n  (B-20)  at  the  numerous  l  values 
required. 

Strictly,  the  leading  product  in  (B-20)  is  a  product  of  N 
principal-value  square  roots.  However,  it  can  be  evaluated 


l)  =  [fT(l-i2Un)  exp  il 


En  +  an) 
1  -  i2lk 


-  «  EZ 

n=l 
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numerically  as  a  single  square  root  of  a  product  of  the  N  terms 
{l-i2£;Xn},  provided  that  the  location  of  this  product  is  tracked 
in  the  complex  plane  from  the  point  1  when  l  *  0.  See,  for 
example,  [21;  page  B-5,  lines  220  -  270], 

If  the  linear  form  in  (B-2)  is  absent,  then  vector  A  =  0, 

A  =  0,  and  an  =  0  for  all  n.  However,  we  still  need  quantities 

en  "  Vn  E  £or  1  1  n  *  N  •  ( B—  23) 

for  the  exponent  in  (B-20),  thereby  necessitating  calculation  of 
eigenvectors  [VnJ.  Only  when  mean  vector  E  is  also  zero  (in 
addition  to  A)  do  just  the  eigenvalues  {XnJ  of  (B-10)  suffice  for 
calculation  of  characteristic  function  f^fS)  in  (B-20). 

In  this  latter  case,  the  matrix  B  C  can  be  considered 
instead,  since  it  has  the  same  eigenvalues  |Xn).  This  follows 
from  the  following  manipulations:  (B-10)  and  (B-ll)  can  be 
written  as 

C  [V1  .  .  .  VN]  =  B- 1  [V1  .  .  .  VN]  diag[Xx  .  .  .  ,  (3-24) 

or 

C  V  -  B-1  V  X  for  1  1  n  £  N  .  (B-25) 

n  n  n 

Therefore 

0  =  (C  -  B"1  Xn)  Vn  =  B~ 1  (B  C  -  Xn  I)  Vn  ( B-26 ) 

or,  since  B  is  positive  definite, 

det ( B  C  -  X  I)  =  0  tor  1  S  n  S  N  .  (B-27) 

'  n 

Thus,  {inl  ar^  the  eigenvalues  of  matrix  B  C  as  well  as  (B-10). 
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CUMULANTS  OF  q 


Before  determining  the  cumulants  of  q,  we  present  a  few 
relations  which  will  enable  easier  manipulations  of  various 
matrices  encountered  below.  From  (B-13),  there  follows 

B  C  =  Q  A  Q- 1  ,  B  C  B  ■  Q  A  QT  , 

B  C  B  C  =  Q  A2  Q"1  ,  BCBCB«QA2QT.  (B-28) 

Therefore 


tr(B  C)  =  tr ( A  Q  Q)  =  tr(A)  =  } _ _  A  , 

n=l  n 


tr | ( B  C ) 2 J  =  tr ( A2  Q-1  Q)  =  tr(Ad 


A2 
n  * 


( B-29 ) 


Now,  if  we  expand  the  natural  logarithm  of  the  general 
characteristic  function  fg(£)  in  (B-20)  in  a  power  series  in  i£, 
we  can  easily  pick  off  the  cumulants  as 

N 

xq(l>  -  Pq  -  ZZ(*„  ♦  ‘l  *  2Vn)  *  tr<B  C)  +  ET  B  E  ♦  2  AT  E  , 

Xq(P)  -  2P_1  (P-1)!  IZ  XP_1  (xn  +  p(cn  +  an)2)  for  p  2  2.  (B-30) 

Here,  we  used  (B-29)  and  (B-22).  In  particular,  the  variance  of 
q  is 


Xq<2>  ■  °\  -  2  EZ  (*2  +  2Xn(En  ♦  a„)2)  - 


2  tr((B  C ) 2 J  +4  (B  E  +  A)T  C  (B  E  +  A)  .  (B-31) 
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This  latter  relation  follows  from  these  manipulations: 

N  2 

IZ  VEn  +  an)  -  (E  +  A)T  a  (E  +  A)  = 

-  (E  +  B-1A)T  Q  A  QT  (E  +  B_1A)  -  (E  +  B_1A)T  B  C  B  (E  +  B_1A)  = 

*  (B  E  +  A)T  C  (B  E  +  A)  .  ( B-32 ) 

Here,  we  used  (B-18)  and  (B-28). 


ALTERNATIVE  DERIVATION  OF  CHARACTERISTIC  FUNCTION  (B-20) 

We  start  again  with  conditions  (B-l)  and  (B-2).  Then,  solve 
(B-10)  for  Q  and  A,  as  before.  Now,  consider  the  linear 
transformation  of  Gaussian  random  vector  X  according  to 

Y  =  QT  (X  -  E)  *  [Yl  .  .  .  yN]T  .  ( B-33 ) 

Then  Nxl  Gaussian  vector  Y  has  mean  Y  *»  0  and  covariance  matrix 


Y  YT  =  QT  (X  -  E ) ( X  -  E)T  Q  -  QT  C  Q  *  A  , 


( B-34 ) 


where  we  used  (B-l)  and  (B-12).  This  diagonal  matrix  means  that 


n 


y  y  =  A  6 
Jn  Jm  " 


n  nm 


( B-35 ) 


that  is,  fyn}  are  zero-mean  uncorrelated  (and  therefore 
independent )  Gaussian  random  variables.  This  is  the  key  to  this 
development. 
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At  the  same  time,  solving  (B-33)  for  X,  we  have,  with  (B-12), 

X  =  E  +  Q-T  Y  *  E  +  B-1  Q  Y  .  (B-36 ) 

Then  quadratic  form  ( B— 2 )  can  be  expressed  as 

q  =  (E  +  B~ 1  Q  Y)T  B  (E  +  B_1  Q  Y)  +  2  AT  (E  +  B_1  Q  Y)  = 

=  YT  Y  +  2  (ET  Q  +  AT  B_1  Q)  Y  +  ET  B  E  +  2  AT  E  = 

rp  fp  fp  fp 

=  Y  Y  +  2  (E  +  A  )  Y+EBE+2AE- 

=  (y2  +  2 ( e  +  a  )y  +  e2  +  2e  a  ]  ,  (B-37) 

L — r  r n  v  n  n'Jn  n  n  n  '  v  ' 

n=l 

where  we  used  (B-12),  (B-18),  and  (B-22).  Therefore,  using 
the  independence  property  derived  in  ( B— 34 )  and  ( B— 35 ) ,  the 
characteristic  function  of  q  is 

_  N  f  dyn  yl 

fa(S)  =  exp(Uq)  =  TT  \  - r  exp  - 

q  n=l  IJ  (2nXn)^  ^  n 

X  exp[u(y2  +  2(  en  +  an)yn  +  *  2e„an)]|  - 

N  f,  \ -*f  (£n  +  °n)2  ?]) 

=  TT  jl1  -  WXn)  expjit  f" T2ir-  -  ilan\\  '  <B-38» 

which  is  equal  to  (B-20). 
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PROPERTIES  OF  EIGENVECTORS  (V  I 

n  ’ 


Consider  the  representation  of  modal  matrix  Q  in  terms  of 
eigenvectors  (Vn)  in  (B-ll).  Then  the  first  relation  in  (B-12) 
and  the  second  relation  in  (B-13)  yield,  respectively. 


T  -1 

V  B  V  =6 
n  m  run 


N 


n=l 


T 

V  V 
n  n 


B 


( B-39 ) 


Similarly,  the  second  relation  in  (B-12)  and  the  fourth  relation 
in  (B-13)  yield,  respectively, 


VT  C  V  -  A  6 
n  m  n  nm 


i  l  Vn  <  -  C_1 
n=l  n 


( B-40 ) 


Thus,  there  are  two  orthogonality  relations  satisfied  by  the 
eigenvectors  (Vn),  namely  the  leading  relations  in  (B-39)  and 
(B-40).  If,  in  addition,  matrix  B  or  C  is  diagonal,  then  the 
trailing  relations  in  (B-39)  and  (B-40)  yield  an  additional 
orthogonality  property. 

In  terms  of  the  components  (B-ll),  the  relations  above  become 


N 

vnk 

k,j-l 

(B  )kj  vmj  **  ^nm 

N 

'  Vnk  ^ 

B)kj  - 

( B-4 1 ) 

and 

1  kj  vmj  ~  ^n  ^nm  ' 

JL 

^ A  vnk  vnj 
n=l  n  J 

o 

i 

*-* 

V_J. 

.  ( B-42 ) 
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GENERAL  SYMMETRIC  MATRIX  B 

All  the  results  above  are  based  on  the  premise  that  matrix  B 
in  the  quadratic  form  in  (B-2)  is  positive  definite.  However, 
there  are  examples  where  this  is  not  the  situation,  in  which  case 
the  corresponding  procedure  for  solution  of  (B-10)  in  appendix  E 
(that  was  mentioned  under  (B-ll))  is  not  applicable.  We  now  give 
an  alternative  procedure  for  obtaining  the  characteristic 
function  of  q  in  (B-2)  for  any  real  matrix  B,  whether  definite  or 
not.  Matrix  B  can  be  taken  symmetric  without  loss  of  generality, 
since  only  the  symmetric  part  of  B  is  active  in  (B-2). 

We  first  observe  that  the  Nxn  covariance  matrix  C  in  (B-l)  is 
always  nonnegative  definite  because 

VT  C  V  =  VT  (X  -  E)(X  -  E)T  V  =  (vT(X  -  E))2  *  0  (B-43) 

for  any  Nxl  real  vector  V.  We  shall  presume  that  C  is  positive 
definite.  Instead  of  solving  (B-10),  we  solve  the  alternative 
generalized  real  characteristic-value  matrix  equation 

B  Q  =  C_1  Q  A  (B-44 ) 

for  Nxn  normalized  modal  matrix  Q  and  diagonal  eigenvalue  matrix 
A.  A  procedure  for  this  solution  is  given  in  appendix  E;  see 
(E-14)  and  sequel. 

Then,  these  solutions  satisfy  [18;  pages  74  -  77] 

QT  C"1  Q  =  I  ,  (B-45) 

QT  B  Q  =  A  s  diag[Aj  .  .  .  ]  .  (B-46) 
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Now  let  linearly  transformed  Nxl  random  vector 

Z  -  Q-1  (X  -  E)  *  [zx  .  .  .  2n]T  .  ( B-47 ) 

Then  mean  Z  *  0  and  covariance  matrix 

z  ZT  »  Q"1  (X  -  E)(X  -  E)T  Q'T  -  Q-1  C  Q'T  *  I  ,  ( B-48 ) 

upon  use  of  (B-l)  and  (B-45).  Solve  (B-47)  for  X,  obtaining 

X  *  Q  Z  +  E  .  (B-49 ) 

Now,  substitute  (B-49)  into  quadratic  and  linear  form  (B-2): 

q  -  (QZ  +  E)Tb  (QZ+E)  +  2  AT  (QZ  +  E)  - 
=  (Z  +  Q~ 1  E)T  QT  B  Q  (Z  +  Q_1  E)  +  2  AT  Q  (Z  +  Q_1  E)  = 

-  (Z  +  B)T  A  (Z  +  E)  +  2  AT  (Z  +  E)  ,  ( B-50 ) 

where  we  used  (B-46)  and  defined  two  auxiliary  N*1  vectors 

E  -  Q~ 1  E  =  [Ej  .  .  .  en]T  ( B-5 1 ) 

and 

A  -  QT  A  s  [ax  .  .  .  Ojj]1  .  ( B-52 ) 

Therefore,  (B-50)  can  now  be  expressed  as 
N 

q  *  )  ( A.  z  +  A  E  +  2  a„)(z^  +  e  )  .  (B-53) 

T  n  n  n  n  n'1  n  n'  '  ' 

n*i 

The  characteristic  function  of  q  can  now  be  readily  found, 
with  the  assistance  of  covariance  property  (B-48),  in  its  most 
compact  form 
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fq(  S )  -  exp(i£q) 


N 


’  (J  dzn  <2">' 

exp[- 

I  »  ,  »1 

exp 

it 


N  A  cj*  +  2e  a  +  i2£cr 
— —  n  n  n  n  n 


n=l 


1  -  i2lK 


n 


( B-54 ) 


In  summary,  the  following  computations  must  be  performed: 
solve  (B-44)  for  Nxn  matrices  Q  and  A;  compute  Nxl  vectors  E  and 
A  by  means  of  ( B—  5 1 )  and  (B-52),  respectively;  and  evaluate  f  It) 
at  desired  £  values  by  use  of  (B-54).  The  only  inverse  matrix 
required  is  that  of  modal  matrix  Q;  the  solution  of  (B-44)  does 
not  actually  require  calculation  of  C-1,  as  will  be  seen  in 
appendix  E,  (E-14)  and  sequel. 

If  mean  vector  E  in  (B-l)  is  zero  and/or  if  linear  form 
vector  A  in  (B-2)  is  zero,  the  corresponding  calculations  in 
(B-51)  and  (B-52)  can  be  circumvented.  If  both  E  and  A  are  zero, 
the  exp  term  in  (B-54)  is  absent  altogether. 

In  this  latter  case,  only  the  eigenvalues  {AnJ  of  (B-44)  need 
be  determined.  If  we  let  Vn  be  the  n-th  eigenvector  (column)  of 
Q,  equation  (B-44)  takes  the  form 

B  V  =  C_1  V  A  for  1  S  n  S  N  .  (B-55) 

n  n  n  '  ' 

This  can  be  manipulated  into 

(B  -  An  C_1)Vn  =  0  ,  ( B-56 ) 

meaning  that 
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det(B  -  Xn  C"1)  -  0  ( B-57 ) 

or,  since  det(C)  t  0,  that 

det(B  C  -  X  I)  «  0  .  ( B-58 ) 

That  is,  {An»  are  the  eigenvalues  of  matrix  B  C.  This  is  a  much 
simpler  numerical  task  than  solving  (B-44)  for  both  Q  and  A. 

If  none  of  the  eigenvalues  |An)  are  zero  or  near  zero,  then 

.  -1  -1  T 

( B-46 )  furnishes  inverse  matrix  Q  »  A  Q  B  and  an  alternative 

to  (B-51),  namely 

-1  T 

E  *  A  1  Q1  B  E  , 

£n  -  j-  B  E  for  1  S  n  S  N  .  ( B-59 ) 

n 

Also,  from  ( B-52 ) ,  regardless  of  the  sizes  of  the  eigenvalues, 

an  *  Vn  A  for  1  S  n  £  N  .  (B-60) 

The  interrelationships  between  Q  and  A  in  (B-44)  and  the 
corresponding  matrices  Q  and  A  in  (B-10)  are  derived  in  appendix 
E,  ( E— 24 )  and  sequel.  The  identity  of  characteristic  functions 
( B-20 )  and  (B-54)  is  also  verified  there. 
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PROPERTIES  OF  EIGENVECTORS  {V  } 

*  n 

Relations  (B-45)  and  (B-46)  yield  the  following  orthogonality 
properties  between  the  eigenvectors: 


T  -1 

V_  C  V  =  6 


V"  B  V  =  k  6 
n  m  n  nm 


(B-61 ) 


Also,  since  (B-45)  and  (B-46)  can  be  alternatively  expressed  as 
C  =  Q  QT  and  B"1  =  Q  A-1  QT,  we  have 


v  v;  =  C  , 

n  n 


C  j-  V  vj  =  B-1 
trrr  Ann 
n=  1  n 


( B-62 ) 


As  a  special  case  of  (B-61),  it  can  be  seen  that  A  =  V  B  V  ; 
therefore,  if  matrix  B  is  nonnegative  definite,  then  An  1  0  for 
all  n.  The  following  example  demonstrates  that  indefinite 
symmetric  R  matrices  can  lead  to  negative  eigenvalues,  even  when 
matrix  C  is  symmetric  and  positive  definite: 


4  -6 


-6  7 


2  3 


( B-63 ) 


the  eigenvalue  and  eigenvector  solutions  to  (B-44)  are 


0  -1 


( B-64 ) 
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CUMULANTS  OF  q 

By  taking  the  logarithm  of  characteristic  function  f  (£)  in 
( B-54 )  and  expanding  in  a  power  series  in  i£,  the  cumulants  of 
quadratic  and  linear  form  q  are  found  to  be 

^  Kt1  +  £n)  +  2 enan  1  £or  m  “  1 

Xq(m)  =  -  > . ( B-65 ) 

(m-l)!  2m-1  XZ  C2[Xn  +  m(XnEn+an)2]  £or  m  2  2. 
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APPENDIX  C-  CHARACTERISTIC  FUNCTION  OF  OUTPUT  y  FOR  GENERAL 
TRANSMITTED  SIGNAL  ENERGIES  {Eki  AND  RECEIVER  WEIGHTS  { AR } 

The  conditional  characteristic  function  of  processor  output 
y,  for  arbitrary  transmitted  signal  energies  {ER)  and  receiver 
weights  {AR},  was  given  in  (27)  as 

K  —  1  k 

fc(S)  *  TT(l  -  U  2ok)  exp | iK  ZZ  "k  rR]  ,  (C-l) 

k— 1  k* 1 


where 


2 

ak 


N  A^ 
No  Ak 


2  ~2 
4  Ak  Ek 


1  "  ^  4  No  Ak 


( C-2 ) 


It  is  important  to  notice  that  the  coefficients  {wR}  in  the 
summation  in  (C-l)  are  complex  and  are  functions  of  argument 
Now,  we  want  to  find  the  unconditional  characteristic 
function  of  y,  after  averaging  over  the  statistics  of  amplitude 
scalings  |rR}.  From  (43)  -  (45), 


2 


M 

zz  x 

m=l 


2 

mk  ' 


( C-3 ) 


where  Gaussian  random  variables 


xmk  -  cm(tk'£k>  +  Wfk*  for  1  S  m  S  H  ,  1  S  k  S  K  .  (C-4 


Therefore,  the  exponential  term  in  (C-l)  becomes 


K 

exp[u  ZZ  wR  rRj  ■  exp[u  XT  W  x] 
k=l 


( C-5 ) 
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where  MKxl  random  Gaussian  real  column  vector 


x  ■  txir"xMi 


•  •  •  Y  •  •  •  Y 

xlk  Mk 


:1K-**XMK] 


(C-6 ) 


and  MKXMK  diagonal  matrix  of  complex  elements 


W  =  diag  [  w^  •  •  •  *w^  •••  ♦ » •  w^  •••  wj^****wK  ] 


( C— 7  ) 


M  terms 


The  MKxMK  covariance  matrix  of  vector  X  is 


Kx  =  Cov { X  j  -  {X  -  X)  (X  -  X)T  ( C— 8 ) 

which  is  known. 

The  unconditional  characteristic  function  of  y  is  available 
from  (C-l)  and  (C-5)  in  the  form 

f  (E)  -  ^TeT  -  rr(l  -  iE  2o2)  exp [ iE  XT  w  x]  - 

k**  1 

-  rr{l  -  iE  2 ol)'1  J  dX  (2„,-mk/2  (det  K .,)■*  X 
k*  X 

X  exp[-  i(X  -  X)T  k~ 1  (X  -  X)  +  H  xT  w  x]  = 

-  rr(l  -  iE  2c2)'1  [det (i  -  iE2  W  Kx))  11  X 
k™  X 

x  exp[i  XT  (I  -  H2  W  Kx)_1  K"1  X  -  \  XT  K”1  x]  .  (C-9) 

The  matrix  W  depends  on  t  through  cu if f icients  {w^J  in  (C-2). 
Therefore,  in  order  to  compute  (C-9),  it  is  necessary  to  invert  a 
complex  MKxMK  matrix  for  each  l  of  interest.  However,  if  mean  X 
is  zero,  only  the  determinant  in  (C-9)  need  be  evaluated. 
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UNCORRELATED  FADING  COMPONENTS  {gro(t,f)} 


In  this  subsection,  we  presume  that  the  fading  components  in 
model  (45)  are  uncorrelated  and  therefore  independent.  This 
means  that,  in  (C-4),  random  variables  xmj  •  •  •  •  'xmK  are  inc*epen- 
dent  of  xni« • • • »xnK  for  m  ^  n.  Reference  to  (C-l)  and  (C-3)  then 
yields  the  characteristic  function  of  output  y  i-n  the  form 


yt>  -  yt>  =Tl[i 


-1 


U  2 ai 


E(S)  , 


(C-10) 


where 


K  o' 

K 

M  2 

E ( l )  =  exp 

U  ^  wk  rk 

k=l  K  KJ 

=  exp 

il  JZZ 

k=l 

wk  C  xmk 
m=i 

M 

TT  exp 
m=l 


K 

u  TZ 

k*l 


w 


k  mk 


(C-ll) 


Now,  let  Kxl  real  random  vector  Xm  and  kxk  complex  diagonal 
matrix  W  be  defined  according  to 

xm  '  >xml  ’  '  '  xmKJT  '  W  =  dia9lwi  •  •  •  »K]  .  (C-12) 


Also,  let  the  covariance  matrix  of  Xm  be  denoted  by  Km,  which  is 
a  kxk  real  matrix,  for  1  S  m  £  M.  Then,  the  m-th  term  in  (C-ll) 
becomes 


exp 


K 

i*  III  w> 

k=l 


X  (det  Km)_J5  exp 


m 


(2n)"K/2  x 


+  itx:  w  x 

m 
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I-H2WK 


Hi  xm  (J  -  “2  w  Km) 


~1  K"1  Xm 
m  m 


1  — T 

4  Xm 

2.  m 


(C-13 ) 


Finally,  the  desired  unconditional  characteristic  function  of 
processor  output  y  follows  from  a  combination  of  (C-10),  (C-ll), 
and  (C-13)  as 


'<£i = ~ u  ^r1  nfKt1  ■ u2  w  kJ)  15 : 

Hi  xm  (J  -  «2  w  Km)'1  C  \  -  I  *m  C  *j}  • 


(C-14 ) 


The  major  computational  burden  here  consists  of  the  inverses 
of  M  complex  Kxr  matrices  for  each  l  value  of  interest.  However, 
if  the  M  covariance  matrices  (K^}  are  all  identical,  the  task  is 
considerably  simplified.  Alternatively,  if  all  the  means  { Xm } 
are  zero,  then  the  evaluation  is  limited  to  M  determinants  of  tne 
complex  Mxm  matrices  I  -  i£2  W  Km-  The  Kxk  complex  matrix  W  is 
defined  by  (C-12)  and  (C-2). 
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APPENDIX  D.  DERIVATION  OF  PROBABILITY  DENSITY  AND 
EXCEEDANCE  DISTRIBUTION  FROM  CHARACTERISTIC  FUNCTION 

The  general  form  in  (136)  and  (158)  of  the  characteristic 
function  of  processor  output  y  can  be  written  as 


We  want  to  determine  the  corresponding  probability  density  func¬ 
tion  and  exceedance  distribution  function.  There  is  a  possible 

numerical  problem  using  Fourier  transforms  on  (D-l)  directly, 

_  k 

since  its  asymptotic  decay  is  only  proportional  to  ^  as  £  -*  °°, 
which  makes  accurate  results  difficult  to  achieve  for  small  K. 

We  will  expand  (D-l)  in  a  power  series  in  ( 1  -  i2£)-*  and 
thereby  obtain  not  only  an  asymptotic  expansion  of  but  in 

fact,  a  convergent  series  which  can  be  transformed  term-by-term 
to  yield  the  probability  density  function  of  y.  Integration  then 
yields  the  exceedance  distribution  function.  Furthermore,  it 
will  turn  out  that  all  the  series  coefficients  are  positive  and 
related  by  recursions  involving  only  positive  terms.  These 
properties  enable  very  accurate  and  efficient  evaluation. 
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We  begin  by  making  the  transformation 


_  1  *  o  t  1  —  2 

z  -  1  -  i2£  '  l2E  ’ - ~ 


(D-3) 


Substitution  in  (D-l)  yields,  after  manipulation. 


fy(S)  =  F  zK  D( z )  E ( z )  , 


where  (using  obvious  abbreviations) 


( D-4  ) 


T  , 
mk 

TP  — 

[Wf(i 

a.  m  1 

mk  1  +  Tmk  '  1 

m,k* 

TmkJ 

D(z)  = 

[Mr  \] 

P1  - 2  r«*)| 

,  E  (  z  ) 

=  exp 

1-2 

2 

(D-5) 


M  gno/j^mk) 

m,k  *  z  rmk 


To  develop  the  series  of  D(z)  in  powers  of  z,  expand 

l  ,M'K 

In  D(z)  =  -  j  TH  ln(l  -  z  r  .)  =  5  !  6  zP  , 


m,  k 


p=l 


( D-6 ) 


where 


i  n 

8  »  s-  >  .  r®\  for  p  2  1  . 

P  2p  t—-r‘  mk 

r  m,  k 


( D-7  ) 


Notice  that  2  0  for  p  2  1  since  Tmk  2  0;  this  latter  property 
follows  from  (D-2)  and  the  fact  that  eigenvalues  of  covariance 
matrices  can  never  be  negative. 

At  the  same  time,  the  expansion  of  In  E(z)  in  (D-5)  is 


l  M  e  * 


p  2P 


(D-8) 
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By  combining  the  results  in  (D-6)  -  (D-8),  there  follows/  after 
some  manipulations. 


ln(D( z )  E ( z ) )  a  zp  ,  ( D-9 ) 

p=0  p 

where 


0  2  r-f  1  + 

m,  k 


mk 


T  ' 
mk 


«  -i&Z 

n  9  L - - 


pP-1 

mk 


P  2  ”'k  'I  +  Tmk>P 


mk 


'mk 


1  +  T 


mk 


for  p  1  1 


(D-10 ) 


Notice  that  Op  2  0  for  p  1  1. 

The  desired  product  for  use  in  (D-4)  is  then  [17;  page  93] 


D(z)  E(z)  *  exp 


«  i  | 

LZ  aD 

p=0  p 


-  ZZ  gp  z1 

p=o  p 


(D-ll) 


where 


g_  =  exp(a  )  ,  g„  =  ^  >  ~  n  a  g„  „  for  p  Z.  1  . 

o  *rx  o'  3p  p  1 — ?  n  3p-n 

n=  l 


( D-12 ) 


It  is  important  to  notice  that  gp  1  0  for  p  i.  0;  that  is,  the 
recursion  in  (D-12)  involves  no  negative  quantities,  thereby 
avoiding  cancellation  error. 

Finally,  the  characteristic  function  in  (D-4)  can  be 
expressed  in  the  desired  power  series  in  z  (=  (1  -  i2£)-*)  as 


U) 


( D-13 ) 
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Now,  we  know  the  following  Fourier  transform  pair  between  a 
characteristic  function  and  a  probability  density  function: 


zn  = 


n-1 


(1  -  i2£ ) n 


for  u  >  0 


2“  (n-1)! 


( D- 14 ) 


This  enables  us  to  write  the  probability  density  function  of  y, 
directly  from  (D-13),  in  the  form 


Py(u) 


uK+p-l  exp(-u/2) 
2K+p  (K+p-1)! 


for  u  >  0  , 


( D-15 ) 


where  scale  factor  F  is  given  by  (D-5).  All  the  terms  in  this 
series  are  nonnegative. 

Integration  on  the  positive  tail  of  density  py  gives  the 
exceedance  distribution  function  of  y: 

00 

oo 

Prob(y  >  u)  =  J  dt  py(t)  =  F  YZ2.  9p  HK+p-l  (f)  for  u  >  0  ,  (D-16) 
u  p-° 

where  F  is  given  by  (D-5)  and 


n  k 

H  (x)  m  exp(-x)  >  fy  for  n  1  0  ,  x  1  0  .  (D-17) 

k=0 


This  latter  sequence  of  functions  is  easily  generated  by  use  of 
the  coupled  recurrences 

Hq(x)  -  exp(-x)  ,  «n(x)  =  Tn(x)  +  Hnl(x)  for  n  M  , 


^Q(x)  *  exp(-x)  ,  Tn(x)  *  Tn_1(x)  ^  for  n  1  1  .  (D-18) 
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The  two  recurrences  derived  above,  namely  (D-12)  and  (D-18), 
utilize  only  nonnegative  quantities.  The  single  negative  term, 
aQ  in  ( D-10 ) ,  is  immediately  converted  to  positive  quantity  gQ  in 
(D-12)  and  aQ  is  not  encountered  again. 

Since  processor  output  y  can  never  be  negative  (see  figure 
1),  we  have,  from  (D-16)  and  (D-17), 

00 

1  *  Prob(y  >  0)  =  F  >  '  g  .  (D-19) 

p=0  p 

This  relation  can  be  used  to  furnish  an  upper  bound  on  the  error 
incurred  by  using  up  to  the  p  =  N  term  in  series  (D-16)  for  the 
exc  idance  distribution  function.  In  particular,  the  error  is 


00 


H» 


(u) 


00 


( D-20 ) 


for  all  u  2.  0.  The  upper  bound  of  1  on  Hn(x)  follows  immediately 
from  definition  (D-17).  Since  scale  factor  F  and  coefficients 
lgpJ  for  0  £  p  S  N  must  be  computed  anyway,  in  order  to  utilize 
(D-16),  error  bound  (D-20)  is  simple  to  compute.  It  says  that 
the  error  at  any  u  is  never  larger  than  the  error  at  the  origin. 

However,  (D-16)  converges  too  slowly  in  some  cases  of  large 
signal-to-noise  ratio;  only  sequence  fg  )  converges  to  zero. 

Xr 

Result  (D-19)  allows  a  modification  to  sum  (D-16)  which  will 
converge  much  more  rapidly.  From  (D-17),  since  Hn(x)  -»  1  as 
n  -*  <»,  we  express  (D-16)  as 
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Prob{y  >  u)  =  F  Zfl  9p  [«K+p_l (5)  -l+l]- 

■  1  -  F^9p[*  -  hk+p-i(i]]  •  <D-21» 

Now,  both  factors  in  the  sum  decay  to  zero  as  p  increases.  In 
fact, 

n  ,  1 

1  "  Hn-l(x)  1  exp( -x )  —  1----  for  n  >  x  .  (D-22) 


This  alternative  in  (D-21)  is  utilized  in  the  program  listed  at 
the  end  of  this  appendix.  Recursions  for  1  -  Hn(x),  which  are 
obvious  modifications  of  (D-18),  are  also  used. 

As  noted  at  the  beginning  of  this  appendix,  form  (D-l) 
encompasses  (136)  and  (158)  in  the  main  text.  It  therefore  also 
covers  (139)  and  (160),  because  these  are  reductions  that  can  be 
obtained,  respectively,  by  setting  all  {emkl  to  zero.  If  all  the 
{emk)  and  {TmkJ  are  zero,  then  all  {ap}  are  zero,  from  (D-10). 
Then,  all  { gp J  are  zero  except  gQ  =  1,  from  (D-12).  In  this  case 
of  noise-only,  series  (D-16)  has  just  the  one  term  HR_^(u/2). 

The  results  above,  for  the  expansions  of  the  characteristic 
function,  probability  density  function,  and  exceedance 
distribution  function,  also  apply  directly  to  forms  (124)  (and 
(126))  as  well  as  (144)  (and  (149))  in  the  main  text,  but  with 
the  following  changes.  For  (124),  define  instead 


T 


n 


e  =  IE  2 
en  N_  En  ' 


F  = 


1  +  Tn) 


( D-23 ) 
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along  with 


KM 


—  _  V — ‘  n 

“o  “  ”  2  2 — =*  1  +  T  ' 
n=  1  ~ 


n 


KM 


a_  *  Tj 

P  2 


tP 

n 


-1 


n=*l  (1  +  T  )! 

n 


n 


n 


P  1  +  T 


n 


for  p  k  1 


The  recurrence  (D-12)  involving  {g^}  is  unchanged. 
For  (144),  define  instead 


2E 


T,  = 


11 


N 


4  ' 


F  = 


ni'  *  'J 


-M/2 


along  with 


1  * 
n  -  i  \ 

o  ~  ~  2  t — r  l  + 


*k 


k=l  *  '  Tk 


a  =  ^ 


,P-1 


k-1  (1  +  T.)P 


M  T  +  -k 

P  k  1  +  T, 


for  p  2.  1 


Again,  the  recurrence  involving  {gp}  is  unchanged. 


(D-24) 


( D-25 ) 


( D-26 ) 
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REFINEMENT 

In  some  cases,  a  refinement  of  the  above  procedure  may  be 
worthwhile,  especially  for  larger  signal-to-noise  ratios.  We 
illustrate  it  with  reference  to  (150),  which  applies  to  M  =  2, 
two  components  in  fading  model  (45).  In  this  case,  we  have 
characteristic  function 


fY(t)  - 


rrf1  -  i2«k) 


-1 


(D-27 ) 


where 


Qk 


1  + 


for  1  i  k  £  K  . 


( D-28 ) 


These  quantities  { Qk }  are  always  larger  than  1  for  nonzero 
signal-to-noise  ratio. 

Instead  of  making  transformation  (D-3),  we  modify  it  to  the 
more  general  form 


Z  =  1  -  i2£Q  ' 


( D-29 ) 


where  Q  can  be  chosen  larger  than  1  if  desired.  Then,  defining 


F 


K 


-1 


TT 

k=l 


for  1  S  k  1  K  , 


( D-30 ) 


characteristic  function  (D-27)  can  be  developed  as 


K  ,  x  1 

-1  V 

00 

fyU)  -  F  zK 

_  K 

=  F  z  exp 

y~2  a  zp 

lp=l  y 

( D-3 1 ) 
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now  with 


QP  " 


K  K 

Eitg  -rz 

k=l  K  k=l 


Qi,  -  Q 


for  p  2.  1 


( D-32 ) 


If  we  take  the  constant  Q  to  have  value 


Q  =  min  Q, 
ISkSK  K 


9 


(D-33 ) 


then  all  the  {t^J  and  [0^}  are  nonnegative,  and  we  still  have 
all  the  desirable  properties  listed  above.  The  value  Q  is  always 
larger  than  1;  it  enables  (D-31)  to  furnish  a  better  fit  to 
(D-27)  with  fewer  terms  than  arbitrarily  forcing  Q  =  1.  Larger 
values  than  Q  can  yield  still  more  rapidly  convergent  series,  but 
at  the  expense  of  some  negative  {g  }.  There  follows,  from  (D-31) 

f  (5)  -  t  zK  fz  gD  ZP  =  P  fZ  9„  (1  '  i2£Q)'K'P  (D-34) 

Y  p=0  P  p=0  p 

for  any  Q;  see  (D-ll)  and  (D-12)  for  the  {9p}«  Expansion  (D-34) 
should  be  compared  with  the  moment  procedure  in  [17;  (102)  and 
(105)],  where  a  different  basis  set,  namely  the  generalized 
Laguerre  functions,  is  used  for  the  series  expansion. 

The  corresponding  exceedance  distribution  function  is 

oo 

Prob(r  >  U)  =  F  n:  9p  HK+p.1(5a)  fo.  u  >  0  .  (D-35) 

Alternatively,  we  have  the  more  rapidly  convergent  form 

OO 

Prob( y  >  u)  -  1  -  F  IZ  9p  [l  -  HK+ p_ j  (jjj)  ]  for  u  >  0  .  (D-36) 
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PROBABILITY  DENSITY  FROM  CHARACTERISTIC  FUNCTION  (163) 

The  characteristic  function  of  normalized  random  variable 
$  =  q/R^0/0)  defined  in  (161)  was  presented  in  (163).  The 
corresponding  probability  density  function  p^(u)  was  then 
evaluated  for  a  few  special  cases  in  (166),  (171),  (175),  and 
(180).  More  generally,  if  we  apply  the  expansion  technique  above 
to  general  result  (163),  we  obtain  the  following  procedure  for 
p^(u):  given  M,  (r^ra^),  Ih^l/  and  threshold  u. 


'  M  -k 

F  -  FT  r(m) 

,m=l 


1  r—  .2 

ao  ~  "  2  ^ — r*  ^m  ' 
m=l 


(D-37) 


-  1  £2  U(m>. . .r_l| -  (rO  -  1  +  p  h*)  for  p  i  1  ,  (D-38) 


r(m) 


gQ  =  exp(aQ)  ,  gp  =  £  i  n  an  gp_n  for  p  M  ,  (D-39) 


h  fu)  =  1  exp( -u/2 ) 

n_ ( n )  -  i 

°  2M/Z  r(M/2) 


hn<U)  "  hp-l(U)  M  -  2  +  ~~2p  f°r  P  2  1  ' 


( D-40 ) 


P+(U)  -  F  9p  hp<U>  for  u  >  0  . 


( D-4 1 ) 


A  function  subprogram  for  this  particular  procedure  is  listed  on 
the  next  page. 
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FUNCTION  SUBPROGRAM  FOR  (D-37)  -  (D-41) 


10  DEF  FNPdf  < DOUBLE  M ,  REAL  U, Rs< * ) , Et asq< * ) ) 

20  DOUBLE  Ms , Ps  !  INTEGERS 

30  ALLOCATE  Q< 1 :  M>  ,  R< 1 : M> , S< 1 : M> , A< 1 : 100) , G<0: 100) 
40  M2*M/'2. 

50  M21*M2-1. 

60  U2*U/2. 

70  F-l. 

80  S*0. 

90  FOR  Ms* 1  TO  M 

100  Rs*Rs<Ms) 

110  T-l./'Rs 

120  E*Etasq<Ms) 

130  GKMs)*l. 

140  R  <  Ms ) *  1 . -T 

150  S  <  Ms ) *E*T 

160  F*F  *Rs 

170  S*S+E 

180  NEXT  Ms 

190  G<0)*EXP(-.5*S> 

200  H*.5*U2''M21*EXP<-U2)/<.SQR<F)*FNGamma<M2>  > 

210  Pdf *G<0)*H 

220  Q-l. 

230  FOR  Ps*l  TO  100 

240  S*0. 

250  FOR  Ms* 1  TO  M 

260  R*R  <  Ms ) 

270  IF  Ps  *  1  THEN  290 

280  CKMs)*Q*Q<Ms)*R 

290  S*S+Q*<R+Ps*S<Ms)> 

300  NEXT  Ms 

310  A<Ps)*.5*S 

320  S*0. 

330  FOR  Ms* 1  TO  Ps 

340  S*S*A<Ms)*G<Ps-Ms) 

350  NEXT  Ms 

360  G<Ps)*G*S/Ps 

370  H*H*U2^<M21+Ps) 

380  T  *G*H 

390  Pdf *Pdf +T 

400  IF  T  <  Pdf  *  1 . E- 1 5  THEN  440 

410  NEXT  Ps 

420  PRINT  "100  TERMS  ARE  INSUF  I C  I  EiJT " 

430  PAUSE 

440  RETURN  Pdf 

450  END 
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A  program  for  the  evaluation  of  (D-21),  by  means  of  (D-10) 
and  recursions  (D-12)  and  (D-18),  is  presented  below.  The 
desired  false  alarm  and  detection  probabilities,  and  P^ 
respectively,  of  the  processor  in  figure  1  are  specified  in  lines 
10  and  20.  Parameter  K  in  line  30  is  the  number  of  signal 
pulses,  while  M  in  line  40  is  the  number  of  components  in  fading 
model  (45).  The  M  power  ratios  ( r ^ ^ }  are  inputted  in  line  50, 
while  the  MK  deterministic  input  signal-to-noise  ratio  measures 
(Dkm/NQ}  are  entered  in  line  60.  The  K  time  locations  JtkJ  and 
frequency  locations  (fkJ  of  the  K  transmitted  signal  pulses  are 
entered  in  lines  70  and  80.  Guesses  at  the  required  average 
random  input  signal-to-noise  ratio  measure  Ej/N0  are  required  in 
lines  90  and  100;  the  program  will  search  for  the  required 
threshold  u  and  required  input  signal-to-noise  ratio  E1/Nq  to 
meet  the  specifications. 

The  normalized  covariance  function  of  the  medium  fading  is 
described  in  function  subprogram  DEF  FNCov(T,F);  currently,  it 
allows  for  exponential  decay  in  time  separation  r  and  frequency 
separation  v,  but  can  be  easily  modified.  It  must  be  noted  that 


covariance  Rmn(t/V)  defined  in  (47)  involves  a  product  of  the 
amplitude- fading  quantities  !gm(t,f)),  not  their  squares.  This 


distinction,  relative  to  covariance  R^^  in  (55)  and  normalized 

2 

covariance  pk-  in  (56)  between  power-fading  variates  {gm(t,f)J, 


must  be  carefully  observed  and  maintained;  furthermore,  it  is 


worthwhile  to  review  (66)  and  (69)  at  this  point,  which  pertain 
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only  to  special  cases  5  and  7. 

Quantities  that  need  be  computed  only  once  and  then  used 
repeatedly  for  different  signal-to-noise  ratios,  such  as  the 
eigenvalues  and  eigenvectors  and  some  auxiliary  arrays,  are 
evaluated  in  the  main  program  and  passed  in  common  to  function 
subprogram  DEF  FNPd  for  the  detection  probability.  Similarly, 
for  speed  and  storage  purposes,  we  have  identified,  for  use  in 
this  subprogram,  the  array  variables  A(p)  =  p  for  p  2.  1  and 
G(p)  =  F  g  for  p  i.  0.  A  separate  listing  for  the  false  alarm 

*r 

probability  in  (188)  is  given  in  function  subprogram  DEF  FNPf . 

The  program  is  listed  in  BASIC  for  the  Hewlett-Packard  9000 
computer.  On  this  particular  device,  the  notation  DOUBLE  denotes 
INTEGER  variables,  not  double  precision.  As  a  numerical  check, 
the  values  printed  out  are  threshold  U  =  38.2583364  and  signal- 
to-noise  ratio  measure  dE  =  10  log10(E1/NQ)  =  13.4744256.  Since 
the  mean  of  processor  output  y  for  noise-only  is  y  =  2K  =  6,  the 
normalized  threshold  is  U/y  =  6.376,  for  this  example  of 
Pp  =  IE-6 . 
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10 

Pf *  1 . E-6 

! 

FALSE  ALARM  PROBABILITY 

20 

Pd*. 9 

j 

DETECTION  PROBABILITY 

30 

K*3 

j 

NUMBER  OF  SIGNAL  PULSES 

40 

M=2 

j 

NUMBER  OF  FADING  COMPONENTS 

50 

DATA  l.,l. 

i 

POWER  RATIOS  r(ro)  FOR 

1 

<  = 

ro 

/  = 

M 

60 

DATA  0,0,0, 0,0,0 

» 

Dkm/No  FOR  1  <=  k  <=  K, 

1 

<  = 

m 

<  = 

M 

70 

DATA  . 1, .3, .5 

i 

TIMES  tk  (SEC)  FOR 

1 

<  = 

k 

<  = 

K 

80 

DATA  .2, .6, .4 

I 

FREQUENCIES  fk  (HZ)  FOR 

1 

<  = 

k 

\  = 

K 

90 

Elr.o0=10. 

I 

E 1 /No  STARTING  VALUE 

100 

E 1  no  1  =  1 . 

i 

El  /■No  INCREMENT 

110 

Ef  *  1 .  E- 15 

i 

TOLERANCE  ON  Pf 

120 

Ed= 1 . E-10 

( 

TOLERANCE  ON  Pd 

130 

PRINT 

140 

PRINT  "Pf  = " ; Pf ; “ 

K  = 

"  ;  K ;  "  M  =  "  ;  M 

150  OPTION  BASE  1 

160  COM  DOUBLE  K,M  !  INTEGERS  (NOT  DOUBLE  PRECISION) 

170  DOUBLE  Ms,Ks,Js  !  INTEGERS 

180  DIM  Rs(5),Dn(10,5>,Ts(19),Fs(10>,Psi (5),Cbar(10, 10) 

190  DIM  U<10, 10) ,V< 10, 10) ,Eig( 10) ,Sq(10,5) 

200  COM  U,Prod(5, 10),Es<5, 10),T(5, 10>,CK5, 10),S<5, 10),B<5, 10) 

210  COM  A( 100) , G<0: 100) 

220  REDIM  Rs< l:M),Dn(lsK, 1 s  M)  , Ts< 1 s K) , Fs< 1 : K) , Ps i < 1 : M> , Cbar < 1 : K,  1 ;  K ) 
230  REDIM  U(liK, 1:K),V(1:K, 1 : K ) , E i g( 1 : K> , Sq( 1 : K , l:M> 

240  REDIM  ProddiH,  l:K),Es(l:M,  1:K),T(1:M,  1  *  K ) 

250  REDIM  Q(1:M, 1 : K) , S< 1 : M, 1 : K) , B( 1 : M, 1 :  K) 

260  READ  Rs(*),Dn(*),Ts(*),Ps(*>  !  Dn(»>  WAS  FILLED  IN  THE  ORDER: 

270  S=0.  !  Dn(l, 1) , Dn( 1 , 2) , Dn( 1 , 3) . Dn(k,m),...,Dn(K,M> 

280  FOR  Ms* 1  TO  M 

290  S=S+Rs(Ms) 

300  NEXT  Ms 

310  FOR  Ms* 1  TO  M 

320  Psi <Ms)=Rs<Ms)/S 

330  NEXT  Ms 

340  FOR  Ks*  1  T0  K 

350  FOR  Js* 1  TO  K 

360  Coo*FNCov(Ts(Ks)-Ts( Js ) , Fs <Ks ) -Fs < Js)  ) 

370  CbarCKs, Js)=Cou  !  NORMALIZED  COVARIANCE  MATRIX 

380  NEXT  Js 

390  NEXT  Ks 

400  MAT  U=Cbar 

410  CALL  Svd(K,K,U(*),V(*),Eig(*)>  !  OUTPUTS:  U<*>, V<*),Eig(*> 

420  PRINT  "EIGENVALUES:" 

430  PRINT  E i g < * ) 
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440  FOR  Ms=  1  TO  M 

450  T=2.*Psi<Ms> 

460  FOR  Ks*  1  TO  K 

4?0  ProdCMs, Ks ) =T*E i gCKs)  !  2  Psi  Eig 

480  NEXT  Ks 

490  NEXT  Ms 

500  FOR  Ks  =  l  TO  K 

510  FOR  Ms= 1  TO  M 

520  Sq<Ks,Ms>=SQRC2.*DnCKs,Ms>> 

530  NEXT  Ms 

540  NEXT  Ks 

550  FOR  Ks  = 1  TO  K 

560  FOR  Ms* 1  TO  M 

570  S=0. 

580  FOR  Js=l  TO  K 

590  S=S+V< Js,Ks>*Sq< Js,Ms> 

600  NEXT  Js 

610  Es ( Ms , Ks ) =S»S  !  emk 

620  NEXT  Ms 

630  NEXT  Ks 

640  Ul  =  l.  !  THRESHOLD  INCREMENT 

650  U0=-LOG<Pf )-Ul  !  THRESHOLD  STARTING  VALUE 

660  CALL  Inversf unct i onl <-Pf , Ef , U0, U1 , U> 

670  PRINT  "THRESHOLD  U  =“;U 

680  CALL  Inversfunc  t  i  on2CPd,  Ed,  Elr>o0,  Elnol ,  Elno> 

690  Db=10.*LGT<Elno> 

700  PRINT  "dB  *" j  Db 

710  END 

720  ! 

730  DEF  FNCov<Tau,Nu>  !  NORMALIZED  COVARIANCE 

740  A=L0G<2.>*.5  !  R1 1 <Tau, Nu)/Rl 1  CO, 0) 

750  B=  i . 3 

760  Cov=EXPC-A*ABS<Tau^-B*ABSCNu)> 

770  RETURN  Cov 

^80  FNEND 

790  ! 

800  DEF  FNPfCU)  !  PROBABILITY  OF  FALSE  ALARM 

810  COM  DOUBLE  K  !  INTEGER 

820  DOUBLE  Ks 

830  U2=U* . 5 

840  Pf*T»EXPC-U2> 

850  FOR  Ks= 1  TO  K-l 

860  T*T*U2/Ks 

870  Pf *Pf  +  T 

880  NEXT  Ks 

890  RETURN  -Pf  !  -  TO  YIELD  INCREASING  FUNCTION 

900  FNEND 

910  ! 
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920  DEF  FNPd(Elno) 

930  COM  DOUBLE  K , M  !  INTEGERS 

940  COM  U, Prod(*> , Es<*> , T(*> , Q(*> , S<*) , BC*} , A<*> , G<*> 

950  DOUBLE  Ms,Ks,Kl,P 

960  To  1  *  1 . E- 1 0  !  RELATIVE  ERROR  OF  SUM 

9?0  FOR  Ms* 1  TO  M 

980  FOR  Ks= 1  TO  K 

990  T<Ms , Ks  >=E lno*Prod  <Ms , Ks ) 

1000  NEXT  Ks 

1010  NEXT  Ms 

1020  F*1 . 

1030  S*0. 

1040  FOR  Ms= 1  TO  M 

1050  FOR  Ks*  1  TO  K 

1060  T=T (Ms , Ks ) 

1070  Q ( Ms , Ks ) =Q= 1 . +T 

1080  S(Ms,Ks)=T/Q 

1090  F  =  F*Q 

1100  S=S+Es(Ms,Ks>/Q 

1110  NEXT  Ks 

1120  NEXT  Ms 

1130  G(0)=EXP(-.5*S>/SQR(F) 

1140  K1*K-1 

1150  U2=U* .  5 

1160  T=EXP(-U2> 

1170  Hl=l.-T 

1180  FOR  Ks* 1  TO  K 1 

1190  T*T*U2^Ks 

1200  H 1 *H 1 -T 

1210  NEXT  Ks 

1220  Pd*l . -G(0)*H1 

1230  FOR  P=1  TO  100 

1240  T*T*U2/(K1+P) 

1250  H1*H1-T 

1260  S*0. 

1270  FOR  Ms= 1  TO  M 

1280  FOR  Ks* 1  TO  K 

1290  Q*Q(Ms, Ks) 

1300  IF  P>1  THEN  1330 

1310  B(Ms,Ks)=B*l . /Q 

1320  GOTO  1340 

1330  B(Ms,Ks>=B=B(Ms,Ks>*S(M3,Ks> 

1340  S  =  S  +  B*(T(Ms,Ks>-'P  +  Es(Ms,Ks>/Q> 

1350  NEXT  Ks 

1360  NEXT  Ms 

1370  fl(P)=. 5*S*P 

1380  S=0. 

1390  FOR  Ms* 1  TO  P 

1400  S*S+A(Ms)*G(P~Ms) 

1410  NEXT  Ms 

1420  G<P)=G=S^P 

1430  De 1 =  G*H1 

1440  Pd=Pd-De 1 

1450  IF  De 1 < Pd*To 1  THEN  1490 

1460  NEXT  P 

1470  PRINT  "100  TERMS  ARE  INSUFFICIENT" 

1480  PAUSE 

1490  PRINT  "Pd  *" ;  Pd,  "  P  =  El/No  =";Elno 

1500  RETURN  Pd 

1510  FNEND 

1520  ! 
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1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

I960 

1970 

1980 

1990 

2000 


SUB  I  river  sf  unc  tionl<Besired,Error,Xl,Del  ,  X2  > 
X2=Xl+De1 
F 1 =FNPf <X 1 ) 

F2*FNPf <X2> 

IF  F2  >=Des i red  THEN  1620 

X 1  =X2 

X2=X2+De1 

F 1  =F2 

GOTO  1560 

IF  F 1 <  =Des i red  THEN  1680 

X2®X1 

Xl*Xl-De1 

F2  =  F  1 

F 1 =FNPf  <X 1 ) 

GOTO  1620 

Xa*Xl 

Xb=X2 

IF  F2-Desired<Desired-Fl  THEN  1770 

T=X  1 

X 1  =X2 

X2  =  T 

T  =  F  1 

F 1  *F2 

F2-T 

IF  HBS<F2-DesiredKError  THEN  1870 
IF  F2=F 1  THEN  1870 

T*<X1*<F2-Desired>-X2*<F1-Desired>)/<F2-Fl > 
T=MRXCT,Xa) 

T=MIN<T,Xb> 

X 1  =X2 
X2*T 
F 1  =F2 

F2=FNPf <X2) 

GOTO  1770 
SUBEND 
! 

SUB  Inversfunct ion2<Desired, Error, XI ,  Del  ,  X2> 

X2=Xl+Del 

Fl=FNPd<Xl ) 

F2=FNPd<X2) 

IF  F2 >  =  Des i red  THEN  1980 
X 1  =X2 
X2=X2+De 1 
F 1  =F2 
GOTO  1920 

IF  F 1 <  =Des i red  THEN  2040 
X2=X  1 
X 1 *X 1 -De 1 
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2010 

F2  =  F  1 

2020 

F 1  =FNPd  <  X 1 ) 

2030 

GOTO  1980 

2040 

Xa=X  1 

2050 

Xb=X2 

2060 

IF  F2-Desi red<De s i red-Fl 

THEN  2130 

2070 

T=X1 

2080 

X 1  =X2 

2090 

X2=T 

2100 

T  =  F1 

2110 

F 1 -F2 

2120 

F2=T 

2130 

IF  ABS(F2-Desi  redXError 

THEN  2230 

2140 

IF  F2=F 1  THEN  2230 

2150 

T=<'X1*(F2-Desired)-X2*(F1 

-Desired))/<F2-Fl) 

2160 

T=MAX(T,Xa) 

2170 

T=MIN(T,Xb) 

2180 

X 1  =X2 

2190 

X2=T 

2200 

F 1  =F2 

2210 

F2=FNPd(X2) 

2220 

GOTO  2130 

2230 

SUBENO 

2240 

! 

2250 

SUB  Svd < DOUBLE  M,  N, REAL  A < * ) , V < * ) , W < * ) ) 

2260 

!  THIS  SUBROUTINE  COMPUTES 

THE  SINGULAR  VALUE  DECOMPOSITION 

2270 

!  OF  AN  ARBITRARY  REAL  MxN 

MATRIX  A:  A  =  U  W  Vt ,  M  >=  N. 

2280 

!  U  IS  MxN ,  V  IS  NxN,  W  IS 

NxNs  W  =  DIAG<D<n>>. 

2290 

ALLOCATE  Rvl(l:N) 

!  NUMERICAL  RECIPES,  PAGES  60-64 

2300 

IF  M >  =  N  THEN  2330 

!  A<*>  IS  OVER-WRITTEN 

2310 

PRINT  "M<N  IS  DISALLOWED" 

2320 

PAUSE 

2330 

DOUBLE  I, J,K,L, Its,Nm, Jj 

!  INTEGERS  (NOT  DOUBLE  PRECISION) 

2340 

G=Sc al e=Anorm=0 . 

2350 

FOR  1=1  TO  N 

2360 

L=  I  + 1 

2370 

Rvl < I >  =  Sc al  e*G 

2380 

G=S=Sc  al e=0 . 

2390 

IF  I >M  THEN  2670 

2400 

FOR  K= I  TO  M 

2410 

Scale  =  Scale  +  ABS<A<K,  I)) 

2420 

NEXT  K 

2430 

IF  Seal e=0 .  THEN  2670 

2440 

FOR  K= I  TO  M 

2450 

Aa=A  <  K , I ) =A  <  K , I)xScal  e 

2460 

S=S+Aa*Aa 

2470 

NEXT  K 

2480 

F=A( 1,1) 

2490 

G=-SQR(S> 

2500 

IF  F<  0 .  THEN  G  =  -G 
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2510  H=F*G-S 

2520  fl(  I ,  I )=F-G 

2530  IF  I =N  THEN  2640 

2540  FOR  J=L  TO  N 

2550  S=0. 

2560  FOR  K  = I  TO  M 

2570  S=S+fl<K, I J) 

2580  NEXT  K 

2590  F*S/H 

2600  FOR  K= I  TO  M 

2610  RCK, J>=fl<K, J)+F*fl<K, I> 

2620  NEXT  K 

2630  NEXT  J 

2640  FOR  K= I  TO  M 

2650  fi<K, I  >*Scal  e 

2660  NEXT  K 

2670  W< I >=Sc  al e*G 

2680  G=S=Sca1 e=0. 

2690  IF  < I >M )  OR  < I =N )  THEN  2990 

2700  FOR  K=L  TO  N 

2710  Scale=Scale+flBS<fld,K)) 

2720  NEXT  K 

2730  IF  Seal e=0.  THEN  2990 

2740  FOR  K=L  TO  N 

2750  fla=fl< I , K ) I , K )/Sc al e 

2760  $-S+fla*fla 

2770  NEXT  K 

2780  F=fl< I , L ) 

2790  G*-SQR<S) 

2800  IF  F< 0 .  THEN  G  =  -G 

2810  H=F*G-S 

2820  R  < I , L  > =F-G 

2830  FOR  K=L  TO  N 

2840  Rvl<K)=fl<I,K)/'H 

2850  NEXT  K 

2860  IF  I =M  THEN  2960 

2870  FOR  J*L  TO  M 

2880  S=0. 

2890  FOR  K=L  TO  N 

2900  S=S+fl<J,K)*fl(I,K> 

2910  NEXT  K 

2920  FOR  K=L  TO  N 

2930  fl< J, K)=fl< J, K)+S*Rvl <K) 

2940  NEXT  K 

2950  NEXT  J 

2960  FOR  K=L  TO  N 

2970  fl<I,K)»fl<I,K)*Scal* 

2980  NEXT  K 

2990  Rnorm=MRX<Rnorm, RBS(W< I ) )+RBS<Rvl < I ) ) ) 

3000  NEXT  I 
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3010  FOR  I =N  TO  1  STEP  -1 

3020  IF  1 >  =  N  THEN  3190 

3030  IF  G=0 .  THEN  3160 

3040  FOR  J=L  TO  N 

3050  V<J, I)=fl<I, J)/fl<I,L)/G 

3060  NEXT  J 

3070  FOR  J*L  TO  N 

3080  S»0. 

3090  FOR  K*L  TO  N 

3100  S*S  +  FKI,K)*V<K,  J> 

3110  NEXT  K 

3120  FOR  K*L  TO  N 

3130  V<K, J>=V<K, J>+S*V<K, I > 

3140  NEXT  K 

3150  NEXT  J 

3160  FOR  J=L  TO  N 

3170  V(I, I>=0. 

3180  NEXT  J 

3190  V<I,I>»1. 

3200  G=Ryl(I) 

3210  L= I 

3220  NEXT  I 

3230  FOR  I=N  TO  1  STEP  -1 

3240  L= I + 1 

3250  G*W<I> 

3260  IF  I >=N  THEN  3300 

3270  FOR  J=L  TO  N 

3280  fi<I,J)*0. 

3290  NEXT  J 

3300  IF  G=0.  THEN  3470 

3310  G*l./G 

3320  IF  I=N  THEN  3430 

3330  FOR  J =L  TO  N 

3340  S=0. 

3350  FOR  K*L  TO  M 

3360  S*S+fKK,  I  J) 

3370  NEXT  K 

3380  F»S/fi<I,I)*G 

3390  FOR  K= I  TO  M 

3400  fl<K,J)=fl<K,J)+F*fi<K,  I) 

3410  NEXT  K 

3420  NEXT  J 

3430  FOR  J*I  TO  M 

3440  fl<J, I)=R<J, I)*G 

3450  NEXT  J 

3460  GOTO  3500 

3470  FOR  J*I  TO  fl 

3480  R<J,I>«0. 

3490  NEXT  J 

3500  fi<  I ,  I  >=fl<  I ,  I  . 
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3510  NEXT  I 

3520  FOR  K  =  N  TO  1  STEP  -1 

3530  FOR  Its=l  TO  30 

3540  FOR  L=K  TO  1  STEP  -1 

3550  Nm=L- 1 

3560  IF  <ABS<Rvl<L>>+Anorm>=f)norm  THEN  3780 

3570  IF  <ABS<W<Nm> >+Anorm>=Anorm  THEN  3590 

3580  NEXT  L 

3590  C=0. 

3600  S=1 . 

3610  FOR  I =L  TO  K 

3620  F=S*Rv 1 < I ) 

3630  Rvl < I )=C*Rvl < I ) 

3640  IF  <ABS<F)+Anorm)=Anorm  THEN  3780 
3650  G=W< I ) 

3660  H=SQR<F*F+G*G> 

3670  W  < I ) =H 

3680  H*1./'H 

3690  C=G*H 

3700  S=-F*H 

3710  FOR  J=1  TO  M 

3720  Ystfl<  J ,  Nm  ) 

3730  Z=A<J,I) 

3740  A< J, Nm>=Y*C+Z*S 

3750  fl< J, I >=-Y*S+Z*C 

3760  NEXT  J 

3770  NEXT  I 

3780  Z*W<IO 

3790  IF  LOK  THEN  3860 

3800  IF  Z>=0.  THEN  3850 

3810  UOO=-Z 

3820  FOR  J=1  TO  N 

3830  V  <  J , K ) =-V  <  J , K ) 

3840  NEXT  J 

3850  GOTO  4390 

3360  IF  1 1 s< 30  THEN  3890 

3870  PRINT  "NO  CONVERGENCE  IN  30  ITERATIONS" 
3880  PAUSE 

3890  X=W<L> 

3900  Nm=K- 1 

3910  Y=W<Nm) 

3920  G  =  Rv 1 C  Nm ) 

3930  H=Rv 1 <  K) 

3940  F  =  <  <Y-Z)*<Y  +  Z)  +  <G-H>*<G-*-H)  )^<2.  *H*Y> 

3950  G=SQR<F*F+1 . ) 

3960  Aa=ABS(G> 

3970  IF  F< 0 ,  THEN  Aa=-Aa 

3980  F=<  <X-Z>*(X  +  Z>+H*<  < Y/CF+Aa) >-H>  >  'X 

3990  C  =  S  =  1 . 

4000  FOR  J*L  TO  Nm 
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4018  I=J+1 

4020  G  =  Rvl < I ) 

4030  Y*W< I ) 

4040  H=S*G 

4050  G*C*G 

4060  Z=SQR<F*F+H*H) 

4070  Rvl < J)*Z 

4080  C*F/Z 

4090  S=H/Z 

4100  F=X*C+G*S 

4110  G=-X*S+G*C 

4120  H  =  Y*S 

4130  Y=Y*C 

4140  FOR  J j  =  1  TO  N 

4150  X=V<Jj,J> 

4160  Z  =  V<Jj,n 

4170  VCJj, J>=X*C+Z»S 

4180  VCJj, I)=-X*S+Z*C 

4190  NEXT  Jj 

4200  Z*SQR<F*F+H*H) 

4210  W<J)=Z 

4220  IF  Z=0.  THEN  4260 

4230  Z-\./Z 

4240  C*F*Z 

4250  S*H*Z 

4260  F=C*G+S*Y 

4270  X=-S*G+C*Y 

4280  FOR  Jj  =  l  TO  M 

4290  Y*fl< J j , J ) 

4380  Z=fi<Jj,I) 

4310  flCJj, J)=Y*C+Z*S 

4320  FKJj, I)=-Y*S+Z*C 

4330  NEXT  Jj 

4340  NEXT  J 

4350  Rvl <L)=0. 

4360  Rvl <K)=F 

4370  W  <  K) *X 

4380  NEXT  Its 

4390  NEXT  K 

4400  SUBEND 
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APPENDIX  E.  SOLUTION  OF  GENERALIZED  EIGENVALUE  PROBLEM 

C  and  B  are  Nxn  real  symmetric  matrices,  and  B  is  positive 
definite.  The  important  alternative  case  where  B  is  not  positive 
definite  will  be  undertaken  in  (E-14)  and  sequel.  We  want  the 
normalized  modal  matrix  Q  and  the  eigenvalue  matrix  A  of  the 
generalized  characteristic-value  problem  [18/  page  74] 
encountered  in  (B-10),  namely 

C  Q  =  B'1  Q  A  .  ( E- 1 ) 

Then,  it  is  known  that  [18;  pages  74  -  77] 

QT  B-1  Q  =  I  and  QT  C  Q  =  A  ,  (E-2) 

which  is  a  simultaneous  reduction  of  two  matrices  to  diagonal 
form.  Alternatively,  when  both  equations  in  (E-2)  are  satisfied, 
then  (E-l)  follows. 

The  first  relation  in  (E-2)  sets  the  scalings  on  the  eigen¬ 
vectors  { Vn ]  of  Q;  in  fact,  it  reads  B-1  Vn  =  1  for  1  £  n  <  N. 

An  alternative  scaling  choice  would  be  to  make  the  eigenvectors 

T 

of  unit  length,  that  is,  Vn  Vn  =  1.  However,  this  would  result 

T  - 1 

in  Q  B  Q  becoming  a  diagonal  matrix  with  non-unity  elements; 

this  latter  alternative  is  not  adopted  here.  Therefore,  the 

eigenvectors  (Vn)  in  (E-2)  do  not  have  unit  length,  that  is, 

T 

vn  vn  /  1  generally.  The  eigenvectors  are  unique  within  their 
polarities . 

A  procedure  for  determining  Q  and  A  is  now  presented.  Solve 
for  the  square  root  matrix  S  in 

B  =  S  ST  ,  (E-3 ) 
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where  S  need  not  be  symmetric.  For  example,  the  lower  (or  upper) 
triangular  square  root  matrix  will  suffice.  Then  there  follows 

-1  -T  -1  T  — 1 

B  =  S  S  and  S  B  S  =  I  .  (E-4) 

Now  compute  Nxn  matrix 

A  =  ST  C  S  ,  ( E-5 ) 

which  is  symmetric.  Solve  the  standard  eigenvalue  problem 

A  U  ■  U  y  (E-6) 

for  normalized  modal  matrix  U  and  diagonal  eigenvalue  matrix  y. 
Then  it  is  known  that  these  solutions  satisfy 

UT  U  =  I  and  UT  A  U  =  y  •  (E-7) 

It  will  now  be  shown  that 

Q  =  S  U  and  A  *  y  (E-8) 

are  the  desired  matrix  solutions  to  (E-l).  Substitution  of 
(E-8)  in  the  first  relation  in  (E-2)  yields 

T-l  TT-l  T 

QBiQ=USBiSU=UIU=I,  (E-9) 

where  we  used  (E-4)  and  (E-7).  Also,  substitution  of  (E-8)  in 
the  second  relation  in  (E-2)  yields 

QTCQ-A=UTSTCSU-y=UTAU-y=0  ,  (E-10 ) 

where  we  used  (E-5)  and  (E-7).  Thus,  both  relations  in  (E-2)  are 
satisfied  by  the  matrix  assignments  in  (E-8). 
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In  summary,  in  order  to  solve  (E-l):  compute  square  root 
matrix  S  in  ( E— 3 ) ;  evaluate  matrix  A  via  (E-5);  solve  (E-6)  for  U 
and  y;  determine  Q  and  A  by  means  of  (E-8).  Notice  there  is 
never  any  need  to  compute  inverse  matrix  B-1;  however,  matrix  B 
must  be  positive  definite  for  the  operation  in  (E-3)  to  succeed. 

We  now  present  an  example  which  illustrates  many  of  the 
quantities  considered  above;  let  matrices 


8 

0 

[2 

1 

B  = 

/ 

C  = 

0 

2 

1 

8 

Then  the  solutions  of  (E-l)  are 


20  01  [2  2 

A  -  ,  Q  =  .  ( E—  1 2  ) 

0  12  J  ll  -1 


The  column  vectors  of  Q  are  not  orthogonal;  in  fact. 


T  ^ 

Q  Q  =  .  ( E— 1 3 ) 

3  5 

Rather,  the  orthogonality  relations  of  (E-2)  may  be  readily 

T  - 1 

verified  to  be  satisfied  m  this  example,  namely  Q  B  Q  =  I  and 
QT  C  Q  =  A. 
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GENERAL  SYMMETRIC  MATRIX  B 


In  attempting  to  solve  ( E- 1 )  for  Q  and  A,  the  square  root  of 
matrix  B  was  required  in  (E-3).  If  B  is  not  positive  definite, 
complex  solutions  would  be  required  there.  Here,  we  will  solve 
the  alternative  generalized  characteristic-value  matrix  equation 
encountered  in  (B-44),  namely 

B  Q  *  C'1  Q  A  ,  ( E-14 ) 

where  covariance  matrix  C  is  real  symmetric  and  positive 
definite,  while  matrix  B  is  real  symmetric  but  can  be  indefinite. 

The  solutions  Q  and  A  will  then  satisfy  the  relations 

QT  C_1  Q  «  I  and  QT  B  Q  =  A  .  (E-15) 


To  this  aim,  first  solve  for  the  square  root  matrix  S  in 


C  =  S  S  . 


( E-16 ) 


Then 


-1  -T  -1  T  -1 

C  =  S  S  and  S  C  S  =  I  . 


(E-17 ) 


Next,  compute  the  Nxn  matrix 

A  =-  ST  B  S  ,  ( E-18 ) 

which  is  symmetric.  Solve  the  standard  characteristic-value 
equation 

A  U  =  U  y  ( E— 1 9 ) 

for  U  and  y,  for  which  we  have  orthogonality  properties 
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UT  U  =  I  and  UT  A  U  =  Y  •  (E-20) 

It  will  now  be  shown  that  the  solutions  to  ( E— 1 4 )  are 

Q  =  S  U  and  A  -  y  •  (E-21) 

Substitution  of  (E-21)  in  the  first  relation  in  (E-15)  yields 
T-l  TT-l  T  T 

Q  C  Q*U  S  C  SU«U  IU«U  U*I  ,  (E-22) 

where  we  used  (E-17)  and  (E-20).  Similarly,  substitution  of 
(E-21)  in  the  second  relation  in  (E-15)  yields 

QTBQ-A=UTSTBSU-Y=UTAU-Y=°  ,  ( E-23 ) 

where  we  used  (E-18)  and  (E-20).  Thus,  both  relations  in  (E-15) 
are  verified  by  solutions  (E-21). 

In  summary,  in  order  to  solve  ( E— 1 4 ) :  compute  square  root 
matrix  S  in  (E-16);  evaluate  matrix  A  according  to  (E-18);  solve 
(E-19)  for  U  and  y,  determine  Q  and  A  by  means  of  (E-21).  This 
procedure  requires  covariance  matrix  C  to  be  positive  definite. 
Inverses  of  matrices  C,  B,  or  S  are  not  required. 

(If  square  root  matrix  S  in  (E-16)  is  taken  lower  triangular, 
then  S_1  is  also  lower  triangular  and  can  be  very  easily  computed 
directly  from  S.  This  enables  calculation  of  inverse  matrix  C_1 
bj  means  of  (E-17),  if  desired.) 
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INTERRELATIONSHIPS  OF  SOLUTIONS 

Here  we  shall  relate  the  solutions  Q  and  A  of  (E-l)  to  the 
solutions  Q  and  A  of  (E-14).  In  particular,  we  maintain  that  the 
specific  relations  are 

Q  -  Q"T  A15  and  A  «  A  .  (E-24) 

To  verify  this  claim,  substitute  (E-24)  into  the  first  relation 
in  ( E-2 )  and  obtain 

QT  B-1  Q  -  I  =  A*5  Q~ 1  B-1  (fT  A*5  -  I  *  A1*  A-1  A**  1  =  0,  (E-25) 

where  we  used  (E-15).  Similarly,  substitution  in  the  second 
relation  in  (E-2)  yields 

QT  C  Q  -  A  =  A**  Q~ 1  C  Q-T  A*5  -  A  =  A*5  I  A**  -  A  «  0  ,  (E-26) 

where  we  also  used  (E-15).  Thus,  both  relations  in  (E-2)  are 
satisfied  by  interrelationships  (E-24). 

The  connection  between  eigenvectors  is  obtained  by  utilizing 
(E-15)  to  give 

Q_T  =  B  Q  A"1  .  (E-27 ) 

Substitution  in  (E-24)  then  yields 

Q  -  B  Q  A'**  ;  that  is,  vn  *  B  Vn  Xn*  *  (E-28) 

The  eigenvalue  relation  is,  directly  from  (E-24), 

Xn  =  An  .  (E-29, 
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Finally,  the  connection  between  vectors  E  and  A  in  (B-18)  and 
the  corresponding  vectors  E  and  A  in  ( B— 5 1 )  and  (B-52)  is  as 
follows:  the  use  of  (E-28)  in  the  first  relation  in  (B-18)  yields 

e  =  V*  E  =  A"**  V*  B  E  =  A*J  e  ,  (E-30) 

n  n  n  n  n  n  '  1 

the  last  relation  based  on  (B-59).  And  substitution  of  (E-28)  in 
the  second  relation  in  (B-18)  yields 

°n  =  Vn  B_1  A  =  K!*  Vn  B  0-1  A  -  A"*5  a  ,  (E-31) 

nn  nn  nn  '  ' 

using  (B-60 ) . 

It  can  now  be  seen  that  the  use  of  (E-29)  -  (E-31)  in 
characteristic  function  (B-20)  immediately  converts  it  to  (B-54). 
This  is  a  direct  confirmation  of  the  fact  that  the  characteristic 
function  of  q  is  unique,  regardless  of  how  derived.  However,  the 
recommended  procedure  is  that  given  in  (B-43)  and  sequel,  because 
it  is  more  general,  allowing  symmetric  matrix  B  to  be  indefinite. 
In  fact,  relation  (E-28)  is  valid  only  when  all  eigenvalues  An 
are  positive.  If  this  is  not  the  case,  although  Q  cannot  be 
found  using  real  arithmetic,  the  technique  summarized  under 
( E-23 ) ,  for  determining  solutions  Q  and  A  of  (E-14),  is  still 
valid  and  can  always  be  used  with  real  arithmetic.  Eigenvalues 
{AnJ  can  take  on  positive  and  negative  values.  The  only 
restriction  is  that  symmetric  matrix  C  be  positive  definite. 
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APPENDIX  F.  STATISTICS  OF  REAL  BILINEAR  FORM 


Suppose  we  encounter  real  bilinear  form  [22;  (13.5-1)] 


b  ■  XT  A  Y  , 


(F-l) 


where  random  vector  X  is  Mxi,  random  vector  Y  is  Nxl,  and 
arbitrary  matrix  A  is  Mxn,  all  quantities  being  real.  We  allow 
random  vectors  X  and  Y  to  be  arbitrarily  correlated  with  each 
other,  but  to  have  zero  means  here;  generalizations  to  nonzero 
means  can  be  made  along  the  lines  of  appendix  B. 

Let  L  =  M+N  and  define  Lxl  vector 


Z  = 


( F-2 ) 


Let  lxl  matrix 


B 


A 

0 


B 


( F-3 ) 


Then  (F-l)  becomes 

b  -  ZT  B  Z  ,  ( F-4  ) 

which  is  a  quadratic  form  in  L  variables;  L  =  M+N.  Matrix  B  is 
symmetric  but  not  positive  definite.  Letting  r  =  min(M,N), 
numerical  investigation  has  revealed  that  B  has  r  positive 
eigenvalues  ;  r  eigenvalues  which  are  their  negatives 

-Pj,...,-pr;  and  the  rest  zero. 

The  mean  of  random  vector  Z  is  zero,  and  its  (arbitrary)  Lxl 
covariance  matrix  is 
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which  is  symmetric.  C  is  also  nonnegative  definite,  because 

VT  C  V  =  VT  Z  ZT  V  =  (vT  z)2  *  0  (F-6) 

for  any  Lxl  vector  V.  We  presume  that  C  is  positive  definite. 

In  order  to  convert  quadratic  form  (F-4)  to  diagonal  form 
with  uncorrelated  variates,  we  first  solve  the  lxl  generalized 
characteristic-value  equation  (see  (E-14)  and  sequel) 

B  Q  =  C"1  Q  A  ( F-7 ) 

for  Lxl  matrices  Q  and  A.  Then,  we  have  properties  [18;  pages 
74  -  77] 

QT  C_1  Q  -  I  and  QT  B  Q  ■  A  ■  diag^  .  .  .  A-L]  .  (F-8) 

Now  let  linearly  transformed  random  vector 

W  >  Q~ 1  Z  »  [wx  .  .  .  wL]  ;  Z  «  QW  .  (F-9) 

Vector  W  also  has  zero  mean.  Its  covariance  matrix  is 

W  WT  -  Q- 1  Z  ZT  Q'T  «  Q- 1  C  O'1  =  I  ,  ( F-10 ) 

upon  use  of  (F-5)  and  (F-8).  Also,  bilinear  form  (F-4)  becomes 

b*ZTBZ=WTQTBQW  =  WTAW  =  £11  '  (F-H) 
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where  we  used  (F-9)  and  (F-8).  Thus,  b  is  a  sum  of  L  weighted 
squares  of  uncorrelated  zero-mean  unit-variance  random  variables. 
Notice  that  not  all  eigenvalues  {X^}  need  be  positive;  some  can 
be  zero  and  some  can  be  negative. 

If  random  vectors  X  and  Y  are  joint-Gaussian,  then  vector  W 
in  (F-9)  is  also  Gaussian,  with  probability  density  function 


exp( -w 


2 

k 


( F-12 ) 


based  on  (F-10).  The  characteristic  function  of  bilinear  form  b 
in  (F-ll )  is  then 


fb(S)  *  exp(Ub)  =  TT(J  dwk  (2n)-Js  exp[-  w£/2  +  UXkw£])  « 

k- 1 


(F-13) 


Since  only  the  eigenvalues  (Xk)  of  A  in  (F-8)  are  required  to 
evaluate  characteristic  function  (F-13),  it  is  not  actually 
recessary  to  determine  the  eigenvectors  Q  in  matrix  equation 
( F-7 ) .  In  fact,  by  a  procedure  identical  to  that  given  earlier 
in  ( B-55 )  -  ( B-58 ) ,  it  may  be  shown  that  {X^J  are  the  eigenvalues 
of  nonsymmetric  Lxl  matrix  B  C.  Reference  to  (F-3)  and  (F-5) 
reveals  that  this  latter  matrix  is  given  by 


T 

xy 


( F— 14) 
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The  rank  of  matrix  B  C  is  generally  R  »  2  min(M,N);  that  is,  the 
product  for  characteristic  function  ffa(£)  in  (F-13)  will  contain 
only  R  terms,  which  is  less  than  L  except  when  M  *  N. 

The  claim  for  the  rank  of  B  C  is  based  upon  the  following 
observation.  Suppose  that  M  £  N;  then  (F-l)  can  be  written  as 


M 


N 


b  - 


m=l  n«l 


m 


a  y 
mn  Jn 


M 

rz 

m»l 


x  v 
m  m 


( F-15 ) 


where  (linearly  transformed)  Gaussian  random  variables 

N 

v  *  )  ;  a_„  y_  for  1  £  m  £  M  .  (F-16) 

m  _ — f  mn  *  n  '  ' 

n*l 

Thus,  b  is  a  sum  of  products  of  just  2M  correlated  random 
variables.  Therefore,  the  characteristic  function  in  (F-13)  can 
have  no  more  than  2M  terms.  Direct  numerical  evaluation  of  the 
eigenvalues  {X^j  of  nonsymmetric  Lxl  matrix  B  C  in  (F-14)  has 
verified  that  its  rank  is  generally  R  «  2  min(M,N),  and  that  the 
R  nonzero  eigenvalues  can  be  positive  and  negative  with  no 
obvious  interrelations.  (If  N  <  M,  the  summation  order  in  (F-15) 
can  be  reversed,  thereby  giving  b  as  a  sum  of  products  of  2N 
correlated  random  variables,  and  hence,  rank  2N  for  matrix  B  C.) 
The  generalization  of  (F-l)  to  the  form 

b  -  XT  An  X  +  XT  A12  Y  +  YT  A21  X  +  YT  A22  Y  ( F-17 ) 

can  be  directly  fit  into  the  framework  of  (F-2)  and  (F-4)  if  we 
generalize  definition  (F-3)  to  LxL  matrix 
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B 


( F-18 ) 


We  can  then  replace  this  B  matrix  by  its  symmetric  part,  since 
only  the  symmetric  part  of  B  is  active  in  quadratic  form  (F-4). 
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APPENDIX  G.  CHARACTERISTIC  FUNCTION  OF  MOST  GENERAL 
COMPLEX  FORM  WITH  FIRST-ORDER  AND  SECOND-ORDER  TERMS 

Let  Z  be  a  complex  Nxl  random  vector  with  real  and  imaginary 
parts  X  and  Y;  that  is,  Z  =  X  +  i  Y.  Let  D^  and  D2  be  complex 
Nxl  constant  vectors.  Let  C^,  C2,  C3,  C4  be  complex  NxN  constant 
matrices,  which  need  not  be  Hermitian  or  symmetric. 

The  most  general  first-order  complex  form  is 

fj  =  dJ  Z  +  D^  Z*  *  (Dj  +  D2)T  X  +  i ( Dj  -  D2)T  Y  —  HT  W  ,  (G-l) 
where 

X 

W  =  .  (G-2 ) 

Y 

H  is  a  complex  2Nxl  constant  vector  and  is  completely  arbitrary; 
that  is,  every  complex  element  of  H  can  be  independently 
specified,  w  is  a  real  2Nxi  random  vector. 

The  most  general  second-order  complex  form  is 

f2  *  ZH  Cx  Z  +  ZH  C2  Z*  +  ZT  C3  Z  +  ZT  C4  Z*  ■  WT  M  W  ,  (G-3 ) 

where 


D1  +  °2 
i(D!  -  d2)'  ' 


Cj  +  ^2  +  ^3  +  ^4  -  C2  +  C3  -  C4) 

.  (G-4 ) 

(Cl  +  c2  -  c3  -  c4)  Ci  -  c2  -  c3  +  c4 

M  is  a  complex  2Nx2N  constant  matrix,  which  need  not  be 
Hermitian.  M  is  completely  arbitrary;  that  is,  every  complex 
element  of  M  can  be  independently  specified. 
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The  pertinent  statistics  of  real  vector  W  are 


,  K  *  Cov(W)  =  (W  -  W)  (W  -  W)T  - 


X  -  X 

f„T  — T 

„T  — T1 

fKxx 

K 

xy 

8 

Y  -  Y 

[Xx  -  xx 

Y  -  Yx]  * 

Kyx 

K 

yyJ 

Here,  w  is  a  real  2Nxl  constant  vector,  while  K  is  a  real  2Nx2N 
symmetric  constant  matrix. 

The  general  complex  form  of  interest  is 

C  *  f2  +  fj  =  WT  M  W  +  HT  W  .  ( G-6 ) 

For  X  and  Y  joint-Gaussian,  W  is  Gaussian  with  probability 
density  function 

p(W)  *  (2n)~N  (det  K)_,s  exp(-  |(W  -  W)T  K-1  (W  -  W))  .  (G-7) 

The  statistical  quantity  of  interest  is,  for  complex  scalar  a, 
the  characteristic  function  of  complex  random  variable  c  in 
(G-6);  in  particular,  the  characteristic  function  is  defined  here 
as  the  average  of  the  exponential 

fc(a)  *  exp(ac )  *  exp|a  WT  M  W  +  a  HT  wj  =  (2ji)  N  (det  K)  **  x 
x  J  dW  exp(-  |(W  -  W)T  K_1  (W  -  W)  +  a  WT  M  w  +  a  HT  w]  = 

=  d_i*  exp(t)  ,  (G-8) 

where 
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d  =  det( I  -  2  a  M  K)  -  det(I  -  2  a  K  M)  , 

1  — T  - 1  —  1  T  - 1 

K  V  +  (I  -  2  a  K  M)  1  k  v  , 

V  =  K"1  W  +  a  H  .  (G-9 ) 

K  is  a  real  2NX2N  matrix,  H  is  a  complex  2Nx2N  matrix  which  need 
not  be  Hermit ian,  W  is  a  real  2Nxl  vector,  H  and  V  are  complex 
2Nxl  vectors,  and  a  is  a  complex  scalar. 

In  general,  we  must  invert  2NX2N  real  symmetric  matrix  K. 
Also,  we  must  invert  2Nx2N  complex  matrix  I  -  2  a  K  M,  which  need 
not  be  Hermit ian,  and  which  depends  on  argument  a.  The  average 
of  the  exponential  in  (G-8)  is  the  type  of  operation  encountered 
in  appendix  C. 

If  X  «  0,  Y  *  0,  Dx  ■  0,  D2  ■  0,  then  W  =  0,  H  -  0,  V  =  0, 
t  =  0,  and  we  need  only  evaluate  complex  determinant 
d  =  det(I  -  2  a  K  M),  which  depends  on  a. 
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SPECIAL  CASE:  Y  =  0 

Then  complex  forms 

fl  "  (D1  +  D2,T  x  s  °T  x  / 

^2  ”  X  (^i  ^  ^2  +  C3  +  C4 )  X  —  X  C  X  , 

C  =  XT  C  X  +  DT  X  .  (G-10 ) 

Matrices  C  and  D  are  complex  and  completely  arbitrary. 

Identify  in  the  subsection  above, 

2N-»N,  W  X  ,  M  ->  C  ,  H-»D,  K  ->  Kxx  ,  (G-ll) 

thereby  getting  characteristic  function 

fc(a)  =  exp(ac)  -  d-*5  exp(t)  ,  (G-12) 

where 

d  =  det{ I  -  2  a  Kxx  C)  , 

1  *  -  \  **  *  +  5  yT  <J  '  2  a  Kxx  C)_1  Kxx  v  ' 

v  =  K~x  x  +  a  D  •  ( G- 13 ) 

Kxx  *s  a  rea^-  NXN  matrix,  C  is  a  complex  Nxn  matrix  which  need 
not  be  Hermitian,  X  is  a  real  Nxl  vector,  D  and  V  are  complex 
Nxl  vectors,  and  a  is  a  complex  scalar. 
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JOINT  CHARACTERISTIC  FUNCTION  OF  REAL  AND  IMAGINARY  PARTS  OF  C 

For  the  general  complex  form  c  in  (G-6),  the  joint 
characteristic  function  of  cr  and  c^  is,  for  real  £  and 

f(S,0  5  exp(  i£  cr  +  iZ,  ci)  .  (G-14) 

But 

5  cr  +  C  ci  =  s(wT  Mr  w  +  Hr  w)  +  C (wT  M£  W  +  H?  w)  = 

=  WT  M  W  +  HT  W  ,  (G-15 ) 

where  real  matrices 

M  s  l  Mr  +  C  Mi  ,  H  h  l  Hr  +  C.  Hi  .  ( G-16 ) 

Therefore,  identifying  a  -»  i,  M  ->  M,  H  -»  H,  in  (G-8),  there 
follows  the  joint  characteristic  function  of  cr  and  c.^  as 

f  ( C ,  C)  *  exp|i  WT  M  W  +  i  HT  w)  =  d_ls  exp(t)  ,  (G-17) 

where 

d  =  det( I  -  i  2  M  K)  *  det(I  -  i  2  K  M)  , 

1  — T  -1—  IT  ~-l 

t  =  -jW  K  W  +  ±V  (I  -  l  2  K  M)  K  V  , 

V  =  K~ 1  W  +  i  H  .  ( G— 18) 
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